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ABSTRACT 

The Fokker-Planck equation has often been used to describe the 
distribution of a particle in a random field. When such a technique is 
applied to a plasma in order to obtain a kinetic equation, certain 
difficulties arise. In the first place, the electric field must be 
determined consistently and the coefficients cannot be found explicitly 
without further assumptions which may not be consistent. Previously 
various attempts have been made to find the coefficients including the 
use of the collision concept by means of the Boltzmann equation and the 
BBGKY theory. In certain cases, this amounts to a manipulation of a 
kinetic equation into the form of a Fokker-Planck equation. Moreover, 
the coefficients are given as averages of functionals of the random 
field. A naive interpretation of these averages, however, gives in- 
correct result. 

In view of these difficulties, a new formulation of the Fokker- 
Planck method is established. It is observed that the coefficients 
actually involve conditional averages of functionals of the random field. 
Furthermore, just as a Fokker-Planck equation can be used to describe 
the lowest distribution, similar equations can also be used to describe 
higher distributions. Such a set of equations is called generalized 
stochastic (or generalized Fokker-Planck) equations following Stratonovich . 
These equations are applicable to systems perturbed by small random 
forces provided sufficient knowledge regarding the unperturbed system is 
available. When the random forces are specified statistically, the 
generalized stochastic equations give immediate results. In this respect, 
of particular interest in plasma application, is the heating of electrons 


IV 


by random electric fields. For the Coulomb potential problem in which 
the electric field has to be determined consistently, it is possible to 
decouple the set of generalized equations by using a cluster expansion. 
The coefficients can then be found explicitly resulting in a kinetic 
equation. The results indicate that the use of a Vlasov equation in 
finding the polarization effect is essentially justified. The polari- 
zation effect is seen to be the consequence of the condition imposed on 
the average in evaluating the coefficient. Unlike the test particle 
theory, however, no test particle is artificially introduced. Finally, 
a similarity is observed between the set of generalized stochastic 
equations and the BBGKY hierarchy. Closer examination reveals that the 
BBGKY hierarchy may be regarded, in a sense, as a special mode of 
generalized stochastic equations. 
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I. INTRODUCTION 

The Fokker-Planck equation has often been used to describe the 
distribution of a particle in a random field. If the distribution is 
spatially homogeneous, the Fokker-Planck equation has the following 
form 1 


3f ( v ,t ) 9 Av f 1_ 9 2 A Av) 2 , f 

at 9v At 2 , 2 At 

av 

where v is the velocity of the particle, Av is the increment of velocity 

in an increment of time At, and f is the velocity distribution. The 

2 

coefficients < -rr > , < — — — > may be found in terms of the random 

field. 1 When applied to the case of a spatially homogeneous plasma, 

certain difficulties arise. In the first place, the electric field 

must be determined consistently and the coefficients cannot be found 

explicitly without further assumptions. Moreover, the polarization 
2 

effect cannot be obtained if the average operator < > is interpreted 

2 

m a naive manner. To find the coefficients Gasiorowicz, et al. , use 
the Vlasov equation for the polarization effect and a Holtsmark 
distribution for the fluctuation of the electric field. It is not 
clear how this can be done since this may lead to inconsistency. 

We observe that the coefficients are actually average values of 
functionals of the random field under certain conditions. Furthermore 
just as a Fokker-Planck equation can be used to describe the distribution 
of the particle, similar equations can also be used to describe higher 
distributions. Such a set of equations is called generalized stochastic 
(or generalized Fokker-Planck) equations. When applied to a plasma, 
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these equations replacing the BBGKY hierarchy can be decoupled by using 
a cluster expansion to determine the coefficients. 

In Chapter II, the generalized stochastic equations for a system 

3 

with only small random forces are derived by modifying Stratonovich' s 
results. For systems having additional forces which are not small, 
reduction is possible through a suitable transformation provided sufficient 
knowledge regarding the unperturbed systems is available. This is con- 
sidered in Chapter III. Some examples in which the random forces are 
assumed to be specified are given in Chapter IV. In Chapter V the 
Coulomb potential problem is considered with generalized stochastic 
equations replacing the BBGKY hierarchy. Using a cluster expansion to 
decouple the equations and freezing the lowest distribution in the 
evaluation of higher distributions, the Fokker-Planck or kinetic equations 
for a plasma are obtained. Although the mathematics involved in this 
case is close to that in the BBGKY theory, the starting point is 
different. It can be seen, however, that the BBGKY hierarchy may be 
regarded as a special form of generalized stochastic equations. 
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II. DERIVATION OF GENERALIZED STOCHASTIC EQUATIONS 
It is known that the fluctuation equation 

x = eF(x,t) (2.1) 

leads to a Fokker-Planck type equation 


3w(x,t) 

3t 


e 


3 _ 

3x 


<F> oj 



0 

K[|^, F ]dt u + 
9x T 


-t 



0 

K[F,F ]dx u + 0(e 3 ) 

T 

-t 


( 2 . 2 ) 


where F = F(x,t), F = F(x,t+x), and K stands for the correlation 

function. The force F as well as the initial value of x are random, and 

3 

to represents the probability density of x at t. Stratonovich gives a 
detailed derivation for a system of equations in the form of (2.1). It 
seems that, however, he has not taken into account of the correlation 
between the force and the initial distribution of x. To account for 
this, we generalize his results and obtain equations similar to (2.2) 
which we shall refer to as the generalized stochastic equations following 
Stratonovich. 

Before we indicate how Stratonovich' s results can be modified and 
generalized, a few remarks concerning notations are in order. In the 
mathematical literature, a random variable is often denoted by a capital 
letter. Thus, X(t) or X^_ represents a family of random variables 
indexed by t , a parameter, whereas lower cases are used to denote real 
numbers. Unfortunately, in the engineering literature, usually no such 
distinctions are made. Thus, the expected value of F(X(0),t) given 
X( 0 ) = x, where F( - ,t) is a random function, can be unambiguously 
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expressed as < F(X(0),t) / X(0) = x> n <F(x,t) / X(0) = x>. It would be 
extremely confusing if not impossible to express such a functional without 
using the capital letters for random variable convention. To compromise, 
we shall follow the engineering literature except in situations where 
it is essential to distinguish a random variable from a real number. 

Thus, to emphasize that x(t) is a random function, we may write X(t) 
instead, and X(t) = x means that the random variable X(t) assumes the 
value of x. 

2.1 Stratonovich 1 s Derivation 

We now review Stratonovich ' s derivation briefly. Consult 

3 

Stratonivch for details. The equation of interest is 


x = eF(x,t) 

with the initial condition 


(2.3) 


x(0) = x r 


(2.4) 


Write the increment by 


x(t) - Xq = eZ^ + e + • . 


= Z 


(2.5) 


and expand F(x,t) as follows 


3 F 

eF(x,t) = eF(x.,t) + e-r— (eZ., + ...) + 
0 8x1 


( 2 . 6 ) 


One then finds 


Z 1 = F ^ x o» t ^ 5 or Z ]/ x o ,1:) = 


, r(x O' t l )dt l 


^2 = fz 1 ' Z 2 (x 0 ,t) 


dt 2M <X 0’ t 2 ) 


F(x 0 ,t 1 )dt 1 etc. 


(2.7) 
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The characteristic function of the random increment 


z(t) = x(t) 


( 2 . 8 ) 


is 


<f> z(t) (u) = 1 + E (lu ? <z n > 


n=l 


nl 


(2.9) 


Inverting 


“z(t) (2/x(0)=x o ) 


1 

2tt 


oo ^ 

-iuz ri _ (iu) „n 

e [1 + E ; — <Z >]du 

. n! 
n=l 


w n 

= [1 + E h <z n >]s< Z ) 


n=l 


( 2 . 10 ) 


where ^ ( z/X( 0 ) -x Q ) is the probability density of Z(t), given X(0) 

- x qi and we shall write as w Z ( t )(z/x 0 ) when confusion is not likely to 


arise. Now 


u x(t) (x/x o } = U) z(t) (x " x o /x o ) 


= Cl + (- ^) <Z n (x 0 ,t)>]6(x-x ) (2.11) 

n=l 


Averaging over x 


0 


X(t) 


(X) = C1 + \ 9x } <2 n (x,t)>]co (x) = [l + L]«o (x) 

n=l (2.12) 


where L is defined by (2.12). Using the more familiar notations 

w(x,t) = u x( t )(x) 

o)(x,0) □ “ X ( 0 )( x ) 
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we have 

u>(x,t) = ( 1+L) ui(x,0) (2.13) 

Taking partial derivative with respect to time in (2.13) and substituting 
o)(x,0) by taking inverse of (2.13), we obtain 


9oj(x ,t ) 
3t 


= L(ltL) 


u(x,t) 


(2.14) 


Expanding L(l+L) , Equation (2.2) is obtained. Details are given by 
3 

Stratonovich. 

2.2 Need for Conditional Averages 

It appears that Stratonovich has assumed that F is not affected by 

the initial specification of X . Otherwise the averages in the charac- 

• -1 

teristic function and hence L(l+L) should be conditional averages. 

This modification is extremely important in the Coulomb potential 
problem as we shall see in Chapter V. In order to show how the modifi- 
cation comes about we first introduce the conditional characteristic 
functions . 

Let Y,Z be random variables. Define the conditional characteristic 

function <j> (u/Y=y) as the conditional average of e luZ given Y = y. 

Zi 

That is 

<}) z (u/Y=y) = <e luZ /Y=y> 

= [e 1U2 to, ( x / Y " y ) d n (2.15) 

J 


= 1 + 7 . 
n=l 


(iu) n 

' n! 


<Z n /Y=y> 


If F does 


not depend on 


explicitly, the results are correct. 


(2.16) 
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where < /Y=y> denotes the expected value given Y = y. Obviously, y 
enters only as a parameter so that the inversion is 


“z (z/Y=y) = 27 


e lu % (u/Y=y)du 


_ 1 _ 

2 7T 


e 1UZ [l + I —j — <Z n /Y=y>]du 
n=l 


= Cl + -4 <Z n /Y=y>]6( z) 

, n ! 3z 
n=l 


(2.17) 


Similarly, if Y' is another random variable, then 

oo Yl 

o^(z/Y=y, Y'=y') = [1 + Z ^j-(- ~) <Z n /Y=y, Y’=y'>]6(z) 


(2.18) 


n=l 


This shows that it is sufficient to replace the averages in the L 
operator in (2.14) with the proper conditional averages. In details, 
suppose we have random variables X, Y, Z such that 

X - Y = Z (2.19) 

Then w x (x/Y=y) = co z (x-y/Y=y) 

00 n 

n [1 + E -4 (- 4-) <Z n /Y=y>]6(x-y) (2.20) 

, n! 8x 
n=l 

Averaging over Y, 

00 XI 

w v (x) = [1 + Z -4 (~ <Z n /Y-x>] a» v (x) 

a -Hi dx I 

n=l 


( 2 . 21 ) 
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If we put X □ X(t), Y = X(0), with to (x) = co(x,t),to (y) = tu(y,0), we have 

A Y 


to(x,t) = [1 + E 
n=l 


l a n _ 

^7 (- ^0 <Z n /X(0)=x>] w(x,0) 


( 2 . 22 ) 


which is Stratonovich 1 s result with the conditional average <Z n /X(0)=x> 
instead of just <Z n >. 

Comparing (2.22) with (2.13), we see that the operator L in (2.13) 
should be modified to L with <Z n /X(0)=x> replacing <Z n >. In order to 
bring out the formal differences due to this modification, let us write 
out the equation corresponding to (2.14) by using the modified operator 
L. We have 

— — = L( 1 + L) to(x,t) (2.23) 

Here 

00 U 

(1 + L) = 1 + E -4 (- <Z n /X(0)=x> (2.24) 

. n: dx 

n=l 

where 

Z = eZ + e 2 Z 2 + . . . (2.25) 

Hence , 

L = - z—~ <Z /X(0)=x> - e 2 <Z /X(0)=x> 

0 X J. dX /. 

2 

+ e 2 (- ■—) <Z 1 Z 1 /X(0)=x> + 0(e 3 ) 


(2.26) 
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and 


Combining 


L( 1 + 


(1+L) 1 h 1-L+ 0(e 2 ) 

= 1 -e(- -||) <Z 1 /X(0)=x> + 0(e 2 ) 

.26) and (2.27) 

_1 = - c ~ <Z 1 /X(0)=x> - e 2 ~ <Z 2 /X(0)=x> 

2 

+ e 2 (- -r— ) <Z, Z /X(0)=x> 

dX 11 

- e 2 (- -£-) <Z /X(0)=x> (- -£-) <Z /X(0)=x> 

3x1 3x 1 

+ 0 (e 3 ) 

= (- 9 |) < eZ 1 + e^' 2 /X(0)=x> 

3 2 2 • 

+ (- ~) <e Z Z /X(0)=x> 

dX 11 

- (- g|)(- g|) <eZ 1 /X(0)=x> <eZ 1 /X(0)=x> 

+ (- -£-) [(- ^-) <eZ /X(0)=x>] <eZ /X(0)=x> 
3x 3x 1 1 

+ 0(e 3 ) 

= e(- -g|o <Z 1 /X(0)=x> 

2 

+ e 2 (- K[Z 1 ,Z 1 /X(0)=x] 


(2.27) 


(2.28) 
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+ e 2 (- ~) <Z 2 /X(0)=x> 

3x 

+ e 2 (- a -) C(- A <Z./X(0)=x>] <Z /X( 0 )=x> 
3x3x1 1 

+0(e^) (2.29) 


where K[A,B/X(0)=x] = <AB/X(0)=,x> - <A/X(0)=x> <B/X(0)=x>. From (2.7) 

we substitute the functions Z , Zg, and find that 


L( 1 + L)" 1 = e( - ^L) <F(x,t)/X(0) = x> 


2 . 3 ' 

+ e ( " 3^ 


K[F(x,t), F(x,t )/X(0) = x]dt. 


2 , 3 

+ e ( 3x 


<F i£^VtI F(x,t )/X(0) = x> dt 
3x1 1 


+ e 2 (- -g|) [(- g|) <F(x,t)/X(0) = x>] 


<F(x,t )/X(0)=x> dtj 


+ 0(e) 


(2.30) 


If F is not affected by the specification of X(0), the third and fourth 
terms on the right hand side of (2.30) combine to give 


2 

e 



K[ 3F(x» t) 


0 


F(x,t 1 )] dt x 


as given by Stratonovich. However, we cannot do so in general. Thus, 
to second order, (2.23) gives the desired equation 
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3io(x,t) 

3t 


-e 


3x 


[<F(x,t)/X(0)=x> 


+ e 


t 



0 


3F(x,t ) 
3x 


F(x,t 1 )/X(0)=x> dt 1 ] 


2 


+ 


2 3 



t 

K[F(x,t), F(x,t 1 )/X(0)=x]dt 1 

0 


+ e 2 -r^- [-T^- <F(x,t)/X(0)=x>] 
3x 3x 


<F(x,t^)/X(0)=x> dt^wtXjt) 


(2.31) 


We remark here that in the cases we consider later, the last term 
contributes to terms of order higher than second and hence can be 
ignored. 


2.3 Generalization to Higher Density Functions and to Higher Dimensions 
Above results can be easily generalized to higher density functions 
and to higher dimensions. 

Let a) , , (x/X' (t 1 ) = x' ) = w(x,t/x' ,t' ) be the probability density 

L ) 

of X(t) = x given X'(t') = x' . Here X'(t') is another random function. 
Equation (2.31) can be applied to the present case, provided we impose 
the condition X'(t') = x' on all averages. This can be seen from 
Equation (2.18), with Y = X(0), Y' = X'(t') where Z is again defined by 
the system (2.3) and 


X(t ) - X(0) = Z 
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Thus , 


w X ( t )(x/X(0) = x Q> X'(t') = x’) 


u z (x-x Q /X(0) □ x Q , X* (t * ) = x' ) 


oo n 

[1 + I (- <Z n /X(0) = x Q , x'(t’)=x>] 6 (x-x Q ) 

n=l 


(2.32) 


Multiplying by ui^ Q ^ (Xq/X 1 ( t ' ) = x') and integrating over x Q , 


a>X(t)( x / x ' ( t ' ) = x ') = uCx.t/x'.f) 


00 a n 

= [1 + l (- ^) <Z n /X(0)=x, X'(t') = x'>] u x(0) (x/X , (t') = x') 

n=l 


" s n 

[1 + l (- <Z n /X(0)=x, X'(t') 

n=l 


= x'>] u (x ,0/x’ ,t' ) 


(2.33) 


Comparing (2.33) with (2.22), we see that the results are formally the 
same provided we impose the condition X'(t') = x' . Thus, (2.33) can be 
rewritten as 

io(X,t/x',t') = (l + L^ f )w(X,0/x',t') (2.34) 

Since the condition X’(t') = x' enters as only a parameter, it does not 
have any effect on the manipulation of L ,(1 + L ,) • Hence, (2.31) 
with the extra condition X'(t') = x' is the equation for ui(x ,t/x' ,t ' ) . 

The results can also be generalized in a straight forward manner 
to multi-dimensional cases. Thus, if instead of (2.3), we have 
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x. = F . ( x ,t ) 
3 3 


j = 1, 2, ..., N (2.3a) 


where 


x = (x lS x 2 , . . . , x N ) 


The generalization of (2.31) is 


0 


9to(x ,t ) 
9t 


e tt— — [<F . /X( 0)=x> + e 
9x . 3 

3 


-t 


9F. 

< F kT /X(0,=X> dt] 

k 


+ e 


2 r 


9x . 9x, 
3 k 


K[F . ,F /X(0)=x]dT 
3 kx 


-t 


+ e 2 <F./X(0)=x>] 

9x. 9x, 3 

3 k 


< F. /X( 0)=x > dx } 03 
kx 


(2.31a) 


Here F_. = F_.(x,t), F^^ = F^^(x,t + i), j ,k = 1,2, ..., N, and the 
convention of summation index has been used. 

To summarize, we have shown in this chapter how Stratonovich' s 
results can be generalized so that the small random forces may not be 
independent of the initial distribution. This allows us to consider 
self-consistent forces other than applied external forces. Furthermore, 
higher density functions also satisfy similar equations. Thus, the 
process can be studied by using such a series of equations which we may 
call generalized stochastic (or Fokker-Planck) equations. They are 
given by (2.31) and the appropriate generalizations. 
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III. GENERALIZED STOCHASTIC EQUATIONS FOR SYSTEMS PERTURBED 
BY RANDOM FORCES 

In the previous chapter, we deal with a system of differential 
equations of the form: 

x = eF 

In other words, the "random force" is small for such systems. In this 
chapter, we wish to study the effect of additional deterministic forces 
which are not small, while assuming we have sufficient knowledge regarding 
the system without the small random forces. 

3.1 General Theorems in Differential Equations 

We begin by collecting some relevant results in differential 
4 5 

equations. * In this section we will use the Einstein summation 
convention. Furthermore, an explicit notation for vectors will not be 
used since the vector property does not play a significant roll except 
in dot products where we use the summation convention. Thus, if we are 
dealing with an N-dimensional space, by x we shall mean the N-tuple 
(x^, ..., x^), and by ^(t,t Q ,x ), we shall mean the N-tuple 
(ij> (t,tg,Xg), ..., ^(t ,t Q ,x) ) . Here t is the scaler time. Thus, 
Equation (2 . 3a) under this convention becomes 

x = eF(x,t) 

which has the same appearance as (2.3). 

Consider the system 

x □ F(x,t) (3.1) 

Under appropriate conditions (such as the Lipschitz condition) on F, 
there is a unique solution x = ij;(t) for each pair (tQ,x Q ) as initial 

Thus, there is a function ^(t,t Q ,x ) such that it satisfies 


conditions . 
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(3.1) and ^(t n ,t Q ,x ) = x Q . Under appropriate conditions (such as 

. 3 F . 

] ] 


(F.~ 3?. 

existence of the first partial derivatives of ^(t ,tg ,x Q ) exist. 


Theorem 1 


Define J = 


3 ( 4 ^, • • • , ^ N ) 

IT" 5 ol 

3(x 1# .... x N ) 


(3.2) 


dJ 8F - 

Then -r— - J -^—1 along each solution and, hence, J = 1 if and only if 
3F. j 

3- 1 = °. 

3x . 

: 

4 

The proof is available m standard literature and will be omitted. 

In the cases we will be interested in, this can also be directly verified. 

Theorem 2 

3F 

In case of an autonomous system (i.e., — = 0), 


<Kt + t, t Q + t, x q ) = ip(t, t Q , x Q ) 


= <P(t - t Q , 0, x Q ) 


(3.3) 


Proof : Let V(t) be a solution through (0,x Q ). Define U(t) = V(t-t Q ). 

Then U(t) is a solution through (t^.x^). But by the definition of 
<Kt,t Q ,x ), we have V(t) = :Jj(t,0,x 0 ), and U(t) = ^( t ,t Q ,x Q ) . Hence, 


i|)(t,t 0 ,x Q ) = U(t ) = V(t-t Q ) □ t-t Q ,0 ,x Q ) (3.4) 


Replacing t by t + t, t Q by t + t in (3.4), we have 


<Kt + t, t Q + t, x Q ) = iKt - t Q , 0, x Q ) 


(3.5) 
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Hence 

'Kt + t, t Q + t, x Q ) = ij»(t - t Q , 0, x Q ) 

= ^’VV 

Theorem 3 

'l'( t » t 0 » x 0 ) = 'Kt,t 1 ,iKt 1 ,t 0 ,x () )) (3.6) 

Proof : ip( t ,t Q ,x Q ) is a solution through ( t^ .^(^ ,t Q ,x ) ) . 

\p( t ,t^ ,t|i( t ,t Q ,Xq ) ) is a solution through the same point, since 

(|»(t ,t .j^Ct ,t 0 ,x Q )) = 1 1 »t ,x Q ) by the definition of \p. By uniqueness, 

these two solutions are identical. 

Theorem 4 

The following is a transform pair for each t, t Q . 

- x = ip(t ,t Q ,y ) 

•< 

. y = iKt 0 ,t,x) 

Proof : Given y, x is uniquely determined by x 

there are y , y^ such that for given x 

x = ^(t,t Q ,y 1 ) 
x = iKt,t Q ,y 2 ) 

Then t|/(t,tg,y^) = i//(t ,t^,y ) since both are solution through (t,x). In 
particular ’K t 0 ,t () ,y 1 ) = y j _ = :/;( t Q ,t Q ,y 2 ) = y ? . By Theorem 3, 

x = \|j(t,t 0 ,y) = 4>(t t 0 ,*(t 0 ,t,x)) 
and the proof is therefore complete. 


(3.7) 

= 4»(t,t Q ,y). Suppose 
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Note that in autonomous case, we may as well choose tg = 0, in 
Theorem 4 and write 

r x = <Kt,0,y) = A(t,y) 

[ y = 1^(0, t,x) = A 1 (t,x) (3.8) 

Corollary 

In case of an autonomous system, then 

A(t,x) = A ^(-tjX) (3.9) 

and 

A(t+o,A ^(t,x)) = A(a,x) (3.10) 

Proof: (3.9) can be readily shown to be true by applying Theorem 2. To 
prove (3.10) we apply Theorem 3. By Theorem 3, 

iKt+a,0,^(0,t,x)) = ^(t+o,t,x) 

= \Ko,0,x) 

which gives (3.10). 

Theorem 3 has the following consequence in partial differentiations. 
Since i|; (t Q ,t+0 ,x) = (t Q ,t ,^(t ,t+o,x) ) we have 


9iJ> 


3x, 


9^ 



t Q ,t t o,x 3x £ 

t Q ,t,^(t,t + a,x) 

3x k 


(3.11) 


t,t + a,x 


In the autonomous case, if A 1 (t,x) is a linear function of x for each 


t, then 


9A . 

3 

9x, 


-1 


-1 

n -1 

9 A . 

9A„ 

] 

£ 

9x. 

9x, 

t + a £ 

t k 


(3.12) 
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where A ^ is defined as in (3.8). Since A ^ is assumed to be linear in 
x, 3 A F /3x is independent of x and hence the subscript x in (3.12) is 
omitted. 

Theorem 5 


3^_.(t 0 ,t,x) 

3t 


+ 


3<K(t ,t ,x) 

F k (x,t) Ix^ ” ° 


(3.13) 


Proof : *.(t 0 ,t,^(t,t 0 ,x 0 )) □ x_. Q by Theorem 3. Hence, 

0 = ^ ^ j (t 0 ,t^(t,t 0 ,x 0 )) 

= al ^(t 0 ,t,4<(t,t 0 ,x 0 )) 

W .(t 0 ,t,^(t,t 0 ,x 0 )) 

+ F k (^(t,t 0 ,x 0 ),t) * 3^- 

Since the left hand side is independent of x Q , we may replace ^(t,t 0 ,x Q ) 
by x (e.g. , by choosing x Q = ^(t Q ,t,x)). (3.13) then follows. More 
generally f( i|>(t ,t ,x) ) satisfies (3.13) where f is arbitrary. Thus 
the partial differential Equation (3.13) with initial condition f(x) has 
f(^(t Q ,t,x)) as its solution. 

Finally we evaluate an expression needed later. Now 


3jKt ,t 
3t 


0 


»y) 


F(i|>(t ,t Q ,y) ,t) 


so that 
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dip 

at 


t,t 0 ,iKt 0 ,t,x) 


= F(iKt,t 0 ,iKt 0 ,t,x)) ,t) 

= F(i|i(t ,t,x) ,t) 

= F(x,t) (3.14) 


3.2 Systems Perturbed by Random Forces 

We now consider the system (3.1) under the influence of small 
random forces. Let the perturbed system be described by 


F. t eF.(x,t) 

: 3 


j = 1 N 


(3.15) 


We shall for simplicity, assume F_. does not depend on t explicitly, 
so that the associated system in the unperturbed state 


x . 
3 


F?(x) 


(3.16) 


is autonomous. 


Let X.(t,t-,x rt ) be the solution of (3.16) with A.(t„,t„,x„) = x„.. 
j 0 0 j 0 0 0 O3 

Since A.(t,t„,x n ) = X .(t-t n ,0,x. ) , we shall choose t„ = 0 and just write 
j 0 0 ] 0 0 0 J 

A^tjX^). Because of Theorem 4, we can define the transformation pair 


r x = t 


t = T 


y . = X . 1 (t ,x) 
3 3 


X. = *.<T,y) 


(3.17) 


where A_. (t,x) = A_.(-t,x) according to (3.19). Equation(3 . 14) is trans- 


formed to 
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dy . 

3 _ 
dx 


3X . 
3 


-1 


at 


dt 

dx 


ax. 

3 


-1 


dx 


3x, 


dt 


k dt_ 
dx 


ax. 1 . ax. 1 dx. 

3 3 k 

at + ax, dt 

k 


ax. 1 n ax, 1 ax. 1 

3 . 3 . _ 3 p 

at k ax, ax, k 

k k 


By Theorem 5 , 


ax. 1 . ax. ' 

-J— + f° — 3- 

dt k 3x, 

k 


0 


so that 

dy. aX. _1 

j— 3- = c -J— F, 
dx 3x, k 

k 


(3.19) 


(3.20) 


Thus system (3.15) goes to the system (3.20) in which the forces are 
small. In (3.20) we should have x and t in terms of y, x using (3.17). 
We shall do this explicitly only for F. Defining 


^■(y,x) = F( X(x ,y ) ,x ) 


(3.21) 


Equation (3. 20) becomes 


dx 



^ k (y»t) 


(3.22) 


Therefore the zeroth order non-random force in (3.15) can be transformed 
away through (3.17). The resulting equation (3.22) has the identical form 
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as the one studied in the last chapter. We shall now apply these 
results to our present system. 

Let tOy(y ) be the density of y at x. From Equation (2 . 31a) : 



(3.23) 

where by ( ) , we mean the function in the parenthesis is to be 

evaluated at a time interval a later. In order to transform (3.23) back 
to variables x and t, we need several relations which we first derive 
in the following. 

Equation (3 . 23) is general and valid even if system (3.16) is non- 
autonomous , provided we take the general transformation (3.7). In such 
case, however, the equation will depend strongly on the initial time. 
Considerable simplification results in a special but important case, 
when (3.16) is linear with constant coefficient so that 

A . "*■( t ,x) = A . ^(t )x . 
i 1 


i = 1, .... N 


(3.24) 
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where the A^ 
designate by ft 


are scalar functions of t, forming a 
■*'(t). Thus in matrix notation 

ft _1 (t) = UT^(t)) 


il 


and (3.17) becomes 

y = ft ^(t)x 

Since A.'V.x) = A_.(-t,x), we can define A „ and 
relations 


A.(x,y) = A„(x)y.. 


and 


x = ft(x)y 


where 


ft( x ) = ( A^j ( x ) ) 


Now we present some interesting properties. They 


A..(x) = A.^(-x) or ft 1 (x) = ft(-x) 
i] !] 




matrix which we shall 

(3.25) 

(3.26) 

ft(x) similarly by 

(3.27) 

(3.28) 

(3.29) 

are 

(3.30) 


(3.31) 
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and 


n(o) = n _ 1 (o) = n(t)n 1 (t) = ft _ 1 (t)ft(t) = i. 


(3.32) 


since A.(0,Xg) = x^ . Here 1^ is the identity matrix. Theorem 1 then 


asserts 


det. ft = 1 = det. ft 


-1 


(3.33) 


3F? 


if the hypothesis = 0 is satisfied. This is always true since the 

X j 

unperturbed system (3.15) is a conservative Hamiltonian system in the 

6 

cases we are interested in. Equation (3 . 12) implies that 


ft \t + o) = ft (t)ft ‘'"(a) 


(3.34) 


= ft _1 (a)ft 1 (t) 


For an arbitrary function H(y,x), we have the following chain rules: 


3H _ 3h (x,t) 3h_(x,t) 

3i " 3t 3 t 3x. (3.35) 

3 


ri X 

3H _ __i 3h(x,t) 

3 7j 9 7j 8x i (3.36) 


where h(x,t) 


a H( A 1 (t,x),t). 


Using (3.13) and (3.31), (3.35) and (3.36) give 


3H 

3t 



F°(x) 

3 


3h 

3x . 
3 


(3.37) 
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9H 



X . . ( t) 
i] 


_9h 

9x . 
1 


(3.38) 


Also, we have 

WyU 1 (t,x),t) = U) x (x,t) (3.39) 

where oo (x,t) is the density of X(t) = x, because of (3.33). Using 
X 

(3.37) through (3.39), equation (3 . 23) can be transformed back to the 
variables x,t. Thus, the left-hand side of (3.23) becomes: 


9w (y ,t ) 9u> (x,t) 0 9u x (x,t) 

~ = 9t~ + F j (x) 9x7 


(3.40) 


The first term on the right side of (3.23): 

-1 

- e 977 < — ^j /Y(0) =y> “y 


- e X, <x7^ F./X(0) = X 1 ( t ,x ) > u (x,t) 

ki9x k i] ] X 


e < F./X(0) = X 1 (t,x)> u. 

9x . 1 


X 


(3.41) 


where we have made use of the fact that X.. is a function of t alone, 

i] 

and that Equation (3 . 32) implies X^X^ = ,X( 0 ) = Y(0). 

For the second term of (3.23): 
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2 _9_ 

3y,- 


da < 


\ in - 1 

9y. 9x, k| 


j k 


/ 


3—^ /Y(0 ) = 


3x^ HI a 


y> 


0 


e 2 \ At)-, 8 


mi 9x 


m 


d0 


= A 1 (t,x)> a) 


Now from (3.21) 


3" AX 1 (t,x),t+a) = F (A(t+a,A 1 (t,x)), t+a) = F (A(a,x),t+a) 


by (3.10). 

Thus (3.42) simplifies to 
0 


.2 _ 9 _ 
9x, 


-t 


3F k -1 

da <^— X iA (a)F^(X(a,x) ,t+a)/X(0) 
i 


X 1 (t,x)> to 


in which we have used (3.32) and (3.34). 
For the third term in (3.23): 


2 9' 


9 9 


0 

'9 A 1 

in' 1 i 

i 

d K 

i v?' 

1 5 

o 

It 

<< 

9x 


-t 

Jv. 

x, / 

J 


= e 2 A ,(t) X .(t) 9 


mi' ' 9x nj 9x J 
m n J 


da K[X. F (t)F k , X^(t+a) F £ (A(a,x 


-t 


= X (t,x)]w. 


t+a)/X(0) 


(3.42) 


(3.43) 


(3.44) 


,t+a)/X(0) 
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2 _a 9_ 

9x. 9x. 
k : 


da K[F k ,X^(a)F £ (X(a,x),t+a)/X(0) 


X (t,x)]a) x ( 3 . 45 ) 


For the last term in (3.23): 
-1 


2 9 


8y i 


3 3A i 

-^ k /Y( 0 ) 


3y. 3x. 

: 


= y> 


° ' 9X - 1 

i 

da < 


3x 


-T 


l I 


3- t o /y(o) = y> » Y 


e 2 X CX <xTJf, /X(0) = X 1 (t,x)>] 

mi3x ni3x lk k 
m n 


da <X_. F (t+a)F^(X(a ,x),t+a) 


-t 


/X(0) = X 1 (t,x)> to 


<F k /X(0) = A_1(t ’ x)> ] 

k 3 


da <X^(a)F 4 (X(a,x),t+a)/X(0) 


-t 


= X 1 (t,x)> to 


(3.46) 


Combining (3.40), (3.41), (3.44), (3.45) and (3.46), 


— + F S ( *- t} asrr = - £ it: <F j /x(0) = A " (t ’ x)> “> 

3 3 


0 


- e 


2 3 


3x, 


3 F 

da k^ F (a)< F^(X(a,x^t+a)/X(0) = X F (t,x)> to^ 

i 


+ e 


2 _9_ 9_ 

9x. 3x. 

k 3 


-1 


- 1 , 


daX^(a) K[F k ,F^(X(a,x) ,t+a)/X(0) = X (t,x)]a» x 


-t 
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+ e 2 < F, / X ( 0 ) 

3x. 3x. k 


X 1 (t,x) >] 


da 


X ]l 


< F^( A(a ,x) ,t+a)/X(0) 



(t,x)> 


(3.47) 


In passing we remark that the basic idea is to transform system 

(3.15) into a system (3.22) with only small forces which will permit us 

to write down a stochastic equation. Such a transformation can be 

effected through the integrals of the unperturbed system, and in our 

case, this amounts to describing the system along its unperturbed orbits. 

As is pointed out in the derivation, Equation (3 . 47) applies in cases 

where (3.16) is linear with constant coefficients and is a conservative 
0 

force. If F is not conservative, there is then need of introducing a 
Jacobian in (3.39) and consequently in (3.47). When (3.16) is non-linear 
but still autonomous, we then have to go back to (3.23). In the most 
general case when (3.16) is non-autonomous the A -transformation in (3.23) 
should be replaced by the more general transformation (3.7). 
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IV. EXAMPLES FOR SPECIFIED RANDOM FORCES 

In this chapter we consider some examples in which the random forces 

are assumed to be specified. Of particular interest in plasma 

application is the heating of electrons by random electric field. 

Sturrock'*'* has considered the heating of electrons without reference to 

.7 

the mechanism for generation of the random electric field. Puri has 
performed an experiment in which a collisionless plasma is heated by 
an electric field of a noisy generator. We obtain similar results by 
using equations developed earlier. The random electric field may be 
due to an incident plane wave and this is considered in the last 
example . 

4.1 Example 1. Charged Particle in a Weak Electric Field 

Consider a charged particle with charge q, mass m, in an electric 

field E(x,t) = - where cf> is the potential. Note that the electric 

oX 

field is not a function of velocity. The equation of motion gives 


x = v 


V 



(4.1) 


If we normalize the equations by' 


x 

V 


t 


* x 0 
^ V 0 


t — 

— v„ 


<t> " ^ <J>Q 


(4.2) 


* We may choose the correlation time of the random field as the 

characteristic time and the thermal velocity as the characteristic 
velocity . 
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where x Q ,v 0 ,cj) 0 are characteristic quantitites not to be confused with 
initial values. Equation (4 . l) in the normalized quantities is then 
given by 


dx 

— = v 
dt 


dv 

dt 


q<J > 0 
( — -) 
v 2 ' 
mv„ 


9^_ 



(4.3) 


q*o 

We express 'weak field" by considering — — = 0(e). Formally then we 

mv 0 

can regard ^-E in (4.1) as 0(e). This corresponds to = v, = 0, 

F = 0, eF 9 a — E . 

1 2 m 


Also 


A,^(t, x, v) = x - vt 


A„ 1 (t, x, v) = 


(4.4) 


or 


fi 1 1 (t) = 


.1 -t. 
■0 1 • 


(4.5) 


Using (3.47), we obtain the stochastic equation for the first example as 


9 to 9(o 

X X 

+ V 


9t 9x 


— | — < E/X(0) = x - vt , V(0) = v > to 
m 9v x 


0 


- (;p) 2 J do(-o) < E(x + va, t+a)/X(0) = x - vt 


m 9v 
V( 0 ) = v > (0 
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. 0 a a 0 

+ (£) 4-— f do(-a) K[E ,E(x h 

m 8v 3x £ 

x - vt, V(0) = v] a) 

+ ( l> 2 Ivl^ j d0 KtE ' E(x + va 

vt, V( 0 ) = v] 

+ (^) 2 |~ [|- < E/X( 0 ) = x - vt 
m 3v 3x 

< E(x + vo, t + o)/X(0) = x - 

+ (^) 2 [|- < E/X( 0 ) = x - vt 

m dV 3v 

< E(x + va), t + a)/X(0) = x 

To compare with Sturrock ’ s results^ we assume the e 
external, with 

< E(x ,t ) > = 0 . 

The electric field is further assumed to be steady 
the statistical sense, i.e. 

< E(x,t) E(x + x', t + a) > 

= < E( 0 , 0 ) E(x ' , o) > 


- vo, t + o)/X( 0) = 

, t + o)/X( 0) = x - 
0 

, V(0) = v >] / do(-o) 

-t 

vt , V ( 0 ) = v > a) 

0 

, V(0) n v >] J do 
-t 

- vt) , V(0) = v > 

(4.6) 

lectric field is 

(4.7) 

and homogeneous in 


= K(x ' , a) 


(4.8) 
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Let S( k ,a)) be the Fourier Transform of the correlation function K(x',o) 
so that 


K(x',o) = If dkdoje^ kx UJ0 '* S(k,m) 

S(k,») = (i^) 2 // dx 'doe -i ^ x ' " “ o) K(x',a) 


Since we have an external field, the conditions on all averages in 
(4. 6) may be dropped. If we are interested in the spatially homogeneous 
case where 


3m 

x 


3x 


0 


(4.10) 


we find Equation (4.6) reduces to 


3co 


- (— ) 2 f do(-a) < E(x + vo , t + a) > to 

m hv J dx 


3t m 3v 


-t 


x 


n 9 a a 

+ (— ) -r— / da < E E(x + vo, t + a) > to 

m 3v 3v J ' 


(4.11) 


-t 


Now 


° 3E 0 

I da(-a) < E(x + vo , t + a) >= J da(-a) < EE(x + va, t+a)> 


0 


8E 

I do(-a) < E — (x + vo , t + a) > 

-{ 3x 


0 


0+ f do < E a — (x + vo, t + a) > 
t 


A 

“ The oj here stands for angular frequency and is not to be confused 
with the probability density function. 
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r 3E 

j do < E (x + va, t + o) > 


0 

J 

-t 


- — f da < EE(x t va, t + a) > 
3v J 


(4.12) 


in which we have used (4.7). The use of (4.12) further reduces the 
stochastic Equation (4 . 11) to 


3(0 


0 


_ - (— ) 2 ~ (w-~ / da < EE(x + va, t + a) >) u> 

3t m 3v 3v i > 


+ (— ) 2 / da < EE(x + va, t + a) > 

m 3v 3v 1 


-t 


or 


3w 0 . 0 s 

Tt" = ( m _) Tv / d ° < EEU + V °’ + 0) > 317 


-t 


3v “x 


( m )2 tv / dt > K( ™- d > 


-t 


3v “x 


Defining 


D( v ,t ) = (— ) 2 / da K(va,a) 

m * 


-t 


(4.13) 


(4.14) 


we obtain finally the desired equation 


3 “x 3 , . 3 

Tt~ = 37 D(v ’ t} 3v W x 


(4.15) 


Here 


D(v, t -*• » ) = ir (% 2 f dk S(k ,vk) (4.16) 

m J 

and by t -*■ <*>, we mean t much larger than the correlation time of E. 


Equation (4 . 15) is a diffusion equation. It has been obtained also by 
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Sturrock. The reader is referred to Sturrock for a discussion of some 
of its properties. 


4.2 Example 2. Charged Particle in a Weak Electric Field with a Constant 
Magnetic Field 

Consider a charged particle under the influence of a constant 
magnetic field perpendicular to its plane of motion. A small random 
external electric field is applied in the transverse direction as 
shown in the figure. For simplicity we assume E is only a function of 
time. The orientation of the magnetic field B and the electric field 
E(t) is depicted in the figure. 

z 



The equations of motion for the present problem are 



v 

x 


, q r 

v = a) v + — E 




0) V 
C X 


+ — E 

m \ 


w 

c 


= S® 


m 


(4.17) 


where 


(4.18) 
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Formally we take E = 0(e) 


(4.19) 


With 



(4.20) 


(4.17) becomes 


*1 = X 3 


x 2 = X 4 


x 3 = “c X 4 + m E x 


X, = -0) x„ + r- E 
4 c 3 m y 


(4.21) 


From the unperturbed system we obtain 


r* 


l (t ,X 10 


X 20 ,X 30’ X 40 ) 


sin a> c t 
X 10 + X 30 u> 

c 


x 40 C c °s a. c t - 1] 
0 ) 

c 


x 3Q [cos u c t - !] X 4 Q sin V 


X 2 (t,X 10’ X 20 ,X 30’ X 40 ) " X 20 + w 


03 


■xC-t-v x x x ) = X„„ cos u t + x sin m t 
X 3^ t,X 10’ X 20 ,X 30’ 40 30 c 40 c 


’ x io ,X 20 ,X 30 ,X 40' 1 


- x 30 Sin V 


x 40 C ° S “c 1 


(4.22) 
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The Corollary (3.9) and Equation (3 . 24) can be used to compute Q 
yielding 


n -1 (t) = 


sin co t [cos co t - 1] 

c c 


[cos co t — 1] sin <o t 

c c 


0 0 cos co t 

c 


0 0 sin to t 

c 


- sm co t 
c 


cos (0 t 
c 


(4.23) 


Substituting into (3.47), 


9co * 3w -x ^ 3u 

— + v • — -r + v x B • t 

9t 9v 


+ e-r 


+ (^r 


3 

9 

9v 

X 

9x 

9 

9 

9v 

X 

3y 

9 

3 

9v 

X 

9v 

X 

9 

_9_ 

9v 

X 

9v 

v 

J_ 

9 

9v 

y 

9x 

9 

9 

9v 

y 

9y 


0 -sm to a 

/ do <j K (o) ( — ) + K (o)[ (cos coo-1)] 

1 I xx (o xy (o c 

-t c J c 


Of sm co a 

/ do i K ( o ) [ - - (cos coo)] + K (o)(- - C -) co 

-t 1 xx “c C xy “c J 


f do < K (o) (cos co o) + K (o) (- sin to a) 
xx c xy c 


— ) 2 t — J da <K (a) (sm <o a) + K (a) (cos co a)l 
m d v 3 v J xx c xv c | 

x v -t J 


-sm co a 

C (a) ( — ) + K [- — (cos co - 1)]> co 

yx co yy co c c | 


yx to c 

J c 


-sm co a 

< (a) ( — H (o 

yy “ c f 
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+ (i ) 2 — J_ 

'm , 9v 9v 

y x 


+ (i ) 2 J-J- 

'm 9v 9v 

y y 


0 

J da 
-t 

0 

| da 
-t 


K (a)[(cos a) a)] + K (a)(-sin a) a ) w 
yx c yy c 


K (a)(sin to a) + K (a)(cos to a) to 
yx c yy c 

(4.24) 


where 


K 

(a) 

= K[E ,E 

xy 

x ya 

K 

(a) 

= K[E ,E 

yx 

y xa 

K 

(a) 

= K[E ,E 

XX 

x xa 

K 

(a) 

= K[E E ] 

yy 

y y° 


(4.25) 


and we have assumed < > - 0 . (4.26) 

(4.24) will be considerably simplified if we assume 


K (a)nK ( a ) = 0 (4.27) 

xy yx 

K (a) = K (a) = K(a) (4.28) 

xx yy 

This is equivalent to saying that the electric field along any 

particular direction has the same statistical property (up to 2nd order). 

Indeed, (4.25), (4.26) and (4.27) imply that < E Q > = 0, < E 0 E 0a > = 

K (a), where E. is the electric field seen in any particular 6-direction, 
xx 6 

If we restrict to homogeneous and isotropic distributions, i.e. 


dui _ 9(jJ 
9x 9y 


0 
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, v . OUJ 

so that "r — = 


to = 

U)(v) 


3to 

V 

X 

3to 

3v 

X 

V 

3v 

3to 

V 

= _x 

3to 

3v 

V 

3v 


/ 2 2 

'v + v 
x y 


(4.24) simplifies to 


a 2 ft 2 n 9 0 

( y + y)(— ) / K(o) COS 0) 0 da I to 

r. 2 „ 2 m i c 

3v 3v -t 

x y l- - 1 


9 3 

2 — — u — — / K(a) cos id o do to 

3u 3u m i c 


where u = ^ mv^ 


Equation (4 . 29) can be rewritten as 


2 2 

_ / 3 ,3 \ ri/ j. \ 

atr = ( — 0 - + - 0 — ) D(t)a) 

9t 3v 2 3 2 

x y 


where 


D(t ) = (— ) 2 f K(c) cos to a do 

m J r 


We note that 


D(t-*») = (3-) / K(o) cos to a da 

m 1 c 

_oo 


= i<a > 2 2 tt s(» c ) 


(4.33) 
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S(u>) 



K(a) cos ojct da 


(4.34) 


Aside from the fact that we allow spatial dependence of the electric 
field in Example 1, the magnetic field is seen to cause a shift in the 
frequency at which the spectral density is to be evaluated. It is clear 
that the solution subject to an initial Maxwellian distribution is 
Maxwellian. In fact if 


w 


-mv 

1 m kT 
2tt k¥ 6 


t 


0 


(4.35) 


then the solution of Equation (4 . 32) is 


-v 


w n 


2 - (T + AT) 
e m 


(4.36) 


2tt -(T + AT) 
m 


where 


AT = — f D(a)da 


(4.37) 


Here k is the Boltzmann constant, T the temperature. The heating rate is 


dAT 

dt 


- m - D(t ) 

K 


(4.38) 


lim dAT 
t -+ °° dt 


2 


= 2tt S(oj ) 
c 


(4.39) 


We have omitted any motion in the z-direction (i.e. along the magnetic 
field). If the correlation between the longitudinal and transverse 
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motions is ignored, the E , E may be evaluated along the rectilinear 

x y 

trajectory in the z-direction. The results agree with Sturrock ' 3 ^ 
Unlike Equation (4 . 6 ), however, Sturrock's diffusion equation is different 
from (4.32) since he uses the guiding center approach. One of the 
disadvantages in using the guiding center approach is that the radius 
of gyration approaches to infinity as the magnetic field decreases 
towards zero. Hence it is difficult to recover the zero magnetic 
field limit from the theory. The present approach does not have such 
difficulties . 

If the magnetic field assumes different values over different 
periods, which are much longer than the correlation time of the electric 

field, we may expect that the result can be readily obtained by 

. . 7 

combining the above results over the intervals. Pqri has considered 

the case with a time varying magnetic field- Except for some minor 

modifications, the approach is similar and straight-forward, so t^at 

we will just briefly indicate how this can be done. Instead of a 

constant B, we let 


a) (t) = a) [1 + A cos (pt + ¥)] 
c c 


where us c , A, p, ¥ are constants. We have 




V 

X 


U (t) 
c 



V 

y 


-oj (t) v + — E 
c x m y 


(4.40) 




(4.41) 
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Instead of the submatrix 

/ 


cos 

U) t 
c 

sin 

0) t 
c 

we have 

AJ 

/ 

VI 

\ 


CO 


s e(t n ,t) -sin e(t ,t) 


>\ 


sin 0(t ,t) cos 0(t n ,t) 


0 : 


/ 


/ \ 


V’ 


(4.42) 


( t t ) 


x 


, j.k = 1,2 


(4.43) 


where 


0(t ,t) □ W c (t-t Q ) + 6 sin(pt Q +'i') - B sin(pt+4') (4.44) 


Aco 


B - 


Note that we do not have an autonomous system, and the more general 
relations (3.7) and (3.11) should be used. The only modification, 
however, is in (3.45) where X.^(o) should be replaced by *F. (t,t+o) so 
that we need only replace cos in (4.29) by cos 0(t,t+o), resulting in 


9tij 

at 


-2 a 2 „ o 0 
r = (-■=■+ -- y)(^) / K(o) cos 6(t,t+o)do 


„ 2 ' 

3v 3v -t 

x y 


(4.45) 
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The result is formally the same with 


0 

D(t ) = J K (a) cos 0(t,t+a)da (4.46) 

-t 


If we merely extend the lower limit -t to -°° in (4.46) D(t) becomes a 
periodic function with a period 2iT/p. Usually the measuring device can 
not respond to such fast changes. Therefore, it seems reasonable to 
take the time average of D(t) for the heating rate. We first expand 
cos 0(t,t+a), using 


cos a sin z = J (a ) + 2 E J a (a) cos 2nz 
u 2n 

n=l 


sin a sin z 


2 E 
n=l 


J 2n+ i(a) sin( 2n+l)z 


where is the Bessel function of the first kind of order n. We find 
that 

cos 0(t,t+a) = cos[a> a + 8 sin(pt t T) - 8 sin(pt + pa + T)] 


cos w a 
c 


^C 1 (pt) C 1 (pt + pa) + C 9 (pt) C 0 (pt + pa) 


2 ^ 2 


sin a^o <C 2 (pt) C^pt + pa) - C 1 (pt) C 2 (pt + pa)f (4.47) 


where 


C 1 (pt) a J Q (8) + 2 ^ J 2n^ ^ COS 2n ^P t + ^ 

n=l 


(4.48) 
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C (pt) = 2 E J 2n+l (e) sin ( 2n+1 )(P t + V) (4.49) 

n=0 

The time average of (4.46) can now be performed by using (4.47). It 
gives 

0 °° 

< D(t) > = (--)^ / K(o) E J? i(6) cos(u +np)o do (4.50) 

t m 1 n c 

_co Yl — — °° 

Here < >^_ means time average. Note that averaging on time is formally 
the same as averaging on 4*. 

Equation (4 . 50) is same as Puri's result. Several points are worth 

notice. Puri uses a randomization period technique. It is not clear 

whether he has taken into account the non-autonomous nature of the 

system. Because of the non-autonomous nature, we find it necessary to 

use the time average technique, which is absent in Puri's treatment. 

Also he has assumed E = 0. In the limit B -► 0 , we see that the 

x 

particle will be accelerated anisotropically . 

Finally, in order to measure the heating rate experimentally, one 
starts out with an ensemble of non-interacting particles with an initial 
distribution. It is pointed out by Puri that if all the particles see 
an identical electric field, the distribution will not be Maxwellian. 

He then argues that if V represents the phase angle of the magnetic 
field along its axis, the particle velocity will consist of a sum of 
' quasi-independent'' modes which enables the use of central limit 
theorem. It is difficult to see how this is done. More basically, if 
only a single realization of the random electric field is carried out, 
it is even more difficult to see the significance. This motivates the 
following example. 
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4.3 Example 3. Plasma Heating by a Plane Wave 

So far we have not considered the origin of the randomness of the 
electric field. We now consider a physical situation that may arise. 

Consider a plane wave propagating in the direction of the magnetic 
field. For simplicity, assume the wave is circularly polarized. The 
equations of motion are taken as 

z = v 

z 

v = 0 
z 

X = V 

X 

• q r 

v = u) v + — E 
x c y m x 



v=-wv+— E (4. 51) 

y c x m y 

The circularly polarized electric field is given by 

E = E A sin(w 0 t - kz) 
x 0 0 

E = E A cos(cj A t - kz) (4.52) 

y 0 0 

where is the frequency and k the wave number of the wave. We have 
ignored the force associated with the magnetic field of the wave. 
Substituting z = z Q + v into the equations for the transverse 
velocities , 


v 

x 


0 ) v + — E.sin((m- - kv A )t - kz ) 
c y m 0 0 zO 0 
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v = -w v + ^ E cos( (to - kv )t - kz.) 
v c x m 0 0 n' 


zO ' 


(4.53) 


Suppose we have an ensemble of non-interacting particles uniformly- 

distributed in space. The initial distribution of longitudinal velocity 

is to be independent of that of transverse velocity. Consider the space 

generated by kz Q = 0, to kz Q = 2tt , for example. (This will roughly be 

a slab if v^t is small.) Let u>(v x ,v ) be the probability density of a 

particle with (v ,v ) in this space. The electric field can be regarded 
x y 

as random with random phase kz Q and random frequency - kv q. For 

any preassigned correlation function K(a), there is a distribution of 

U) A - kv or a distribution S(v .) of v . so that the electric field 
0 zO zO zO 

Q 

is a stationary process with this correlation function. In fact 


<E > = <E > = 0 
x y 

< E (t) E (t + o) > = < E (t) E (t + o) > 

xx y y 


= 7 E o J cos[(w 0 ~kv z0 )o]S(v z0 ) dv z0 

_oo 

= K(o) 


(4.54) 


(4.55) 


Here S is seen to be proportional to the spectral density of E (t). 
The cross correlations are not zero but antisymmetric , 

< E (t) E (t + o) > 
x y 



45 


= - < E (t) E (t + a) > 

y * 


= | E l j sinC^ - kv zQ )a S(v zQ ) dv zQ 

—OO 


Applying (4.24), we have the following diffusion equation. 


(4.56) 



and the heating rate is 


2 

m< 




(4.59) 


Note that the maximum occurs when ^c + ^0 equals the most probable 

k 

velocity along the direction of propagation. 

4 . 4 Summary 

In this Chapter, we have applied our previous results to several 
situations, in which a charged particle experiences a random electric 
field. By considering a noninteracting ensemble of charged particles, 
it is found that heating is possible with a random electric field, 
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and the heating rate is in general proportional to the spectral density 
of the random field. In general, the presence of a magnetic field has 
the effect of a shift in the frequency at which the spectral density 
is to be evaluated. This is particularly clear in the last example, 
in which the origin of the randomness of the field is also considered. 
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V. GENERALIZED STOCHASTIC EQUATIONS APPLIED 
TO THE COULOMB POTENTIAL PROBLEM 

In the introduction we have already pointed out that the application 
of the conventional Fokker-Planck equation to the Coulomb potential 
problem involves some difficulties. For one reason the coefficients can- 
not be determined explicitly without the use of some model or other theory. 
From a theoretical point of view, this may lead to inconsistency. On 

the other hand, the kinetic equation of a plasma based on the BBGKY theory 

. . 9 10 

can be manipulated into the form of a Fokker-Planck equation. ’ Aside 

from the fact that we wish to derive the kinetic equation using the 

Fokker-Planck equation approach, this does not give a general exact form 

of the coefficients in terms of the random force. 

Early attempts in finding the coefficients introduces the collision 

10 12 2 

concept by means of the Boltzmann equation. * Gasiorowizc et al. , 

use the Vlasov equation and a Holtzmark distribution. Following 

13 

Bogoliubov's method, Tolmachev introduces a chain of linked distributions 

at different times for the Coulomb potential problem. The results are 

14 

divergent at extreme distances of interaction. Temko makes use of a 
ternary correlation to insure convergence at large distances. His results 

9 

have been further refined by Tchen. These methods are primarily based 
on the Liouville equation or the BBGKY theory, and are especially 
formulated for the Coulomb potential problem. As a result, these methods 
may not be suitable for systems with more general forces. 

In this chapter gneralized stochastic equations are used for the 
lowest and higher distributions. With a cluster expansion, the equations 
are decoupled and the coefficients are determined. 
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Kinetic Equation 


The dotted line indicates some relation between the generalized stochastic 
equations and Lionville equation. This point is discussed in the final 
chapter. 

5.1 Formulation of the Coulomb Potential Problem 

We consider a system of N identical electrons in a volume V, with 
a uniform background of immobile positive charge. Let be the spatial 
coordinate of the jth electron, v_. the velocity and E_. the electric 
field seen by the electron. The equations of motion are then 


V. = M- eE.) 
3 m 3 


j - 1 , •••> N 


(5.1) 
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where -e and m are the charge and mass of the electron. Introducing the 


potential 



(5.2) 


we have 



d<P . 

dx. 

J 


(5.3) 


Normalizing Equation (5.1) by 


x . 

1 


= x .r 
-] 


0 


v . 

1 


v.v 

-1 


0 


m 


t = 



Equations (4 . D become 


(5.4) 


dx . 

-=1 = v. 
dt ^ 



2 

dt mv Q 


3<j> . 

rJ-) 

9x . 
“3 


(5.5) 


Here we take r^ an effective range of the force, v^ as the thermal 
velocity, with 
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(5.6) 


4> . - E —x 3 

] k/jlv ? k 


(5.7) 


The characteristic time is — and should be such that the higher 

V 0 

distributions relax in a time of a few characteristic times, enabling us 
to drop the initial conditions and extend the limit of integration t to 
00 as in the later development. We assign 


4b 


(5.8) 


mv. 


to be a small number, 
the thermal energy. 
Let 


This means the interaction is weak compared to 


£ = (x,v) 


(5.9) 


be a six dimensional vector. Let D _ _ (£ ,t ; ... ; £.,t.) be 

h x ? ... , y_. xi ] 3 

the probability density of (x (t^j'v (t^)) = E,^, ... , (x^ (t.),"v^ (t.)) 


jki'sft - 1 - - 1 - . . - 1 

□ 1,2, ... . Here y maps the set {1, ... , j) into (1, . 


y 3 3 
, N}. 


We consider only symmetrical distributions, i.e., if M maps (1, ... , N} 
onto (1, ... , N} in a 1-1 manner, then 


°M(p 1 ), ... , M(y.) ( Vi ; » 5 j’ t j ) 


V •• 


("I, ,t. ; ... ; "t . ,t . ) 


, u, l’ 1 


3 3 


(5.10) 


For a plasma, r_ is usually taken as the Debye length and the 
characteristic time is taken as the reciprocal of the plasma frequency. 
kit 

To keep our notations simple and comparable to those in the literature, 
we will not use capital letters as in Chapter I. 
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Since the system (5.1) is invariant under the transformation M, this 
. . . . ^ 

is equivalent to requiring D 1 N (? l’°’ *** » ? N’°^ t0 be s y mmetr>ical - 

In view of the symmetry, we will suppress the subscripts y^, ... y. 


whenever they are distinct. Thus 


D(C 1 ,t 1 ; ... ; Sj.tj) - D Xj ^ j^l^l* ’*• » 


D M(i), ... , M(j) u i’V ••• ; 5 j j ■ 


(5.11) 


Let 


V 






. ;Cj » t j/^ 1 »tqi 




V * 


»V y l 


2^9 ••• 9hj < ( ^l»"^l » •*' 9 5 j 9^ j 9 j .. 


' , ... •** *H* t k ) 




(5.12) 


be the conditional probability density of 


( \ (t l , ’\ (t l ,) = «!• • 


. ,(x (t.),v (t.)) = i. 

y j 3 y j 3 3 


given 






where 


j = 1,2, ... and k = 1,2, ... 
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Again, when the subscripts y^ . . . , y^,y^, ... ,y^ are distinct, we will 

just write as D(^ 1 »t 1 ; ;5j ... ;C^,t^.). 

Since we have a large volume V, it is convenient (particularly in 

3 D ( t ) 

the case of spatial homogeneity, i.e., — = 0) to write 

3x 

D(t,t) = i f^t.t) (5.13) 


More generally, we let 


D(f 1 ,t; . 


itj.t) = <±P f^. 


• ,t) , j < N 


(5.14) 


Similarly let 


D( » t ^ ; 


5t j^i ,1: i’ 


rf' t' ) 

’^k’ k ; 


( V ) j g j/k ( '^l ,t l’ 


;Z- 


^,t'), j + k < N (5.15) 


16 


We now introduce the cluster expansion. We shall write down only 
the first few relations explicitly. Define P, T, Q as follows: 


D(t 1 ,t 1 i5 2 ,t 2 ) = D(E 1 ,t 1 )D(5 2 ,t 2 ) + P(t 1 ,t 1 ;? 2 ,t 2 ) (5.16) 


3 ^ 3 ^ 

D( "l’ t l ; ^2’ t 2’^3’ t 3 ) = D^.t.) + V D(5 i ,t.)P(C j ,t.;T k ,t k ) 

j <k 


i=l 


+ T ^i» t i’^2 ,t 2’^3 ,t 3^ 


(5.17) 
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D(C 1 ,t 1 ;C 2 ,t 2 ;? 3 > t 3’ ? 4» t 4 ) = 11 D(5 i> t i ) 

i=l 


-i. 


+ V D(C.,t.)D(5 j ,t.)P(5 kj t k ; V)l ) 
i<j ,k< l 




i<j ,k<& 


t Q^i» t i’^2 ,t 2 ’^3 ,t 3’^4 jt 4^ 

(5.18) 

where E' means the summation indices i, j, k, i are to be all distinct. 

In the later development, we require P to be of first order, T and 
Q of higher order than P. To zeroth order, since the trajectories of 
the particles are independent, this implies that 


D(5-^5T^;5 2 »"t2’^3’^'3’^4 9 ^'4^ 


ir D(f.,t i ) 
i=l 


(5.19) 


and 


1 = 1 


(5.20) 
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Similarly, we have 


g j/k^l ,t l’ * 


;? j ,t j /5 l ,t l’ ’ 


• = yl 11 


i=l 


1 l 


tt f ( E , . ,t . ) j+k < 4 
...111 J 
1 = 1 


(5.21) 


to the zeroth order. We will not employ special notations for 
D 

It should be noted in such cases, D 
In fact, to zeroth order. 


(£ ,t ; ... t.) when the p.'s are not all distinct. 

,...,Pj±i HD ) 

may contain 6-functions. 

] 


4-^ ••• ,4. 


D -t- t •~F t •'?" • "?' t' +■ ' -■?' 'i 

1,2,3,4,1,2,3,4^1’ L 1’^2’ 2’^3’ 3’S’ 4 5 ^1’ l’^2’ 2’^3’ 3’S’ 4 ; 


: TT D( 5 . ,t . )6 (x ! -"x .-v. ( t '. -t . ) )6( v! - V - . ) 
11 11111 11 


1 4 4 

(y) tt fj_(£p >t^ )6(x|-lL -v\ (t|-t^ ) )6("v! -‘v. ) 
i=l a ' J ' 


(5.22) 


where 


£• = ,v. ) 

l i i 


c: = <*.*) 


We are now ready to write down some stochastic equations. Let us 
define the following quantities. 


eE . = eE . (x . ,t . ) 

:a □ ot i n 


(5.23) 
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N 9<j>(x . ,x. (t . )) 
y d k p 

3x. 

k^p Da 


• , N 


(5.24) 


4>(x) = 

Ixl 


(5.25) 


E . □ E. 0 (x.+v.a,t.+a) 

pBa pB p p ’ p 


(5.26) 


5.(t) = (x.(t),'v.(t)) 
D D 2 


j = 1,2, ... , N 


(5.27) 


K. = (x. $.) 
D D P 


j = 1,2, ... , N 


(5.28) 


rT. = (x.-v.t.,v.) 
D DID] 


j =1,2, ... , N 


(5.29) 


where the subscripts a, prefer to the ath and gth component of a three 
dimensional vector. 

Taking j □ 1 in system (5.1), we can write down the stochastic 
equation for D('t 1 ,t ) or f^"^,^). This is formally the same as in 
example 1 (note that however the potentials are defined somewhat 
differently). To allow for a three dimensional space, refer to Equation 
(2.31a). The result is: 


3 f l (C l’ t l ) 3f l e /?ro)='n>f 

3t + V la 9x m 3v. la 1 11 

1 la 1 
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2 ^ 
e 9 


2 9v, 
m la 


9E 

do( _ 0 )<_^ E i fj<j / ti( o ) = fr, > f, 

lp 


2 

e 

9 

9 

2 

9v_ 

9 v _ . 

m 

la 

16 l 

2 , 

e 

9 

9 

2 

9 v. 

9x, . 

m 

la 

13 J 

2 

e 

9 

[ 3 

2 

9 v_ 

l 9x, „ 

m 

la 

16 

2 

e 

9 

[ 3 

2 

9 

L 9v 1 0 

m 

la 

13 


dcKtE 1 a- I W !: 1 (0) = "l> f l 


0daK[Ei a ,Ei ea /;i(°) = q 


do(-a) <E leo /? 1 <0) = V f- 




d “ <E 1 Bc /{ 1 <0) = V f l 


(5.30) 

where a, 3 are to be summed from 1 through 3, and we have replaced the 
lower limit -t, of integration by -<*>, since we are interested in time 
long compared to the correlation time. 

According to the results in Chapter I, D( 5-^ »t^/5 2 >*2 > a 

similar Equation as D(£^ t^) so that gj_/q( ’^2^ sat i s f i- es (5.30) 

provided the extra condition 


C 2 ( V = 5 2 


(5.31) 
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is imposed on all averages , so that , 

^i/i^i’V ^’V + v 9 g i/i 


9t, 


la 9x 


la 


m 9v <E la / ^l (0) “ " ^2* g l/l 

la 


e 2 9 


2 9v n 
m la 


9E, 


d<J (. 0) <_, E 1Bo / {l (0) = n 1>{2 (t 2 ) 
IP 


? 2 > g l/l 


2 

e 

9 

9 

2 

3 v_ 

3v_ _ 

m 

la 

IB J _ 

2 

e 

9 

9 

2 

9 v. 

9x. . 

m 

la 

16 

2 

e 

9 

r 9 

2 

3 v_ 


in 

la 

16 

2 

e 

9 

r 8 

2 

3 v. 

L 9v, „ 

m 

la 

16 


d ° KCE la’ E lga /C l (0) = V«2 (t 2 } = 5 2 ]g l/l 


da(-a)K[E la ,E lgo /5 1 (0) = '^ 1 .5 2 (t 2 ) = ^ S l/1 


3, /?.( 
la 1 


2 s 2 


da(_a)<E 16a /C l (0) 


h 1 ,S 2 (t 2 ) 


la 1 1 ”2 2 ' 


? 2 > g l/l 


^2 >] 


do<E 13a /? l (0) 


n 1 ’^ 2 ( t 2 ) ? 2 > g l/l 


(5.32) 


Thus, if we know the various conditional averages, we can find f and 
g l/l‘ ^ ac1: That These conditional averages of E^ are functions 

of f^, and perhaps higher density functions and must be determined 
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consistently. This is clear from Equation (5 . 3) which defines E . Indeed, 
if we make use of the definition of the various density functions, we 
find that 


< E. /^(O) 
la 1 


n l > 


< Z ■ 

j=2 


N 3<t»(x 1 -x.(t 1 )) 


3x 


la 


/sfO = V 


M 3$(x -X.) g 

I, 8x U j . t 1 /<t l’ Q)de j 

j =2 la 


(N-l) 


3*|>CSc 1 --ft 2 ) ^ ^ 

g l/l ( WV° )d5 2 

la 


(5.33) 


where we have made use of the assumption that the particles are 
indistinguishable. Similarly, 


' E i« E ieo /? i <0) = V 


N 3<*»(x -t ) 3<f>(x +v .o-x' ) g . , ^ j. 

J. —kr 1 H^-^VvW°'v 0)d ^ 

la lp V 

k,j>l 


N 


9 cj> ( x -it . ) 3<|)(* 1 +^ 1 o-«') . * J . l 

E - 1 ^ 


k=j*l la 


IB 


(N-l)(N-2) 


3<}> (■$-$) 3 4> ( x +v a-x' ) 


3x 


1 2 


la 


13 
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+ ( N— 1 ) 


a <f>(x’ 1 -x 2 ) 


3ifi(x 1 +v 1 a 


-*;) 


3x 


la 


3x 


IB 


" D 2,2/l (5 2, t i ; 


=2 >t ,+a/rv^ jOM^dCg 


(5.34) 


Using (5.21), (5.22), (5.33) and (5.34), we can estimate the orders 
of magnitude of the various terms in Equation (5 . 30). The ratios of these 
terms are found to be 

1st 
term 

1 

(5.35) 


2nd 


3rd 


4th 


5th through 8th 


1 : n o r ? ( V mv o 2) 


n o r o ( V mv o 2)2 


(n o r o ) 2 ( V mv o 2)2 


Here we have taken N, V to be large, with a finite ratio 

N _ 

V " n 0 (5.36) 

3 

We will consider two cases. In the first case, n^r^ 'u 1, so that 
(5.35) is of the form 

1 : 1 : e : e 2 : e 2 (5.37) 

2 

Keeping terms up to e , we will show that (5.30) results in the Fokker- 

3 1 

Planck equation. In the second case, we consider the limit n Q r Q ->■ — , 
which leads us to the kinetic equation. In both cases, we assume 
spatial homogeneity, i.e. , 

8f l ( ^l’ t l ) 


0 
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5.2 The Fokker-Planck Equation 
In this case we consider 

n^r^ ^ 1 (5.38) 

and keep terms of (5.30) to second order. According to (5.33), (5.34) 

and (5.37), we need g^^ to first order in (5.33), gj/^ as we ll as 

D to zeroth order in (5.34). The last two can be obtained from 

2 , 2/1 

(5.21) and (5.22). Thus our primary concern is to find g^^ to first 
order. Note that the ordering in (5.37) applies to (5.32) as well. Thus 
to first order, gj_/j_ satisfies 


3g l/l ( ^l ,t l /? 2 ,t 2 ) 

at. 


8g 


t v 


1/1 


la 9v 


la 


£ 

m 3v 


3_ 

la 


E i. /E i (0) 


n 1 ,5 2 ( t 2 ) ^2 > g i/l 


(5.39) 


By symmetry of the density functions as well as that of system (5.1) 
(under the transformation M) , 


9g l/l,^ ^ ^ , 3g l/l _ e 

(? 2 ,t 2 /5 l ,t l ) + V< la3v 
2 2a 


m 3v 


2a 




- n 2 ,5 1 (t 1 ) = 5 X > g 1 / 1 


(5.40) 


The result is expected since it is simply Equation (5 . 32) or Equation 
(5. 39) for the second electron. By first setting t^ = 0, then 
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9g l/l ( '^2’ t 2 / ^ ,0) 


+ V 2oi 3v g l/l m 9v <E 2a /5 2 (0) = n 2’ 5 l (0) 
2 2a 


= V g l/l 


.a -a 


Now <E 2a^2^°^ - r, 2 , ^l^^ = n l > ds Eoun ^ in a way similar to (5.33). 


e <E 2a /5 2 (0) = n 2 .e 1 C0) = v 


3<Kx -x. ) g , 

} 9x7 — (?j.t 2 /n 1 ,0;n 2 ,o)dc j 

• l ]=3 2a J J 


9(j)(x -x|) ^ ^ 

+ - 9x 2a D;l / 1> 2 ( ^l’ t 2 /n i’ 0;n 2 ,0)d ^l 


f 3 4 >(x 2 -x 3 ) ^ ^ ^ 

n 0 9x7 g l/2 (? 3 ,t 2 /n l’° ;n 2’ 0)d ^3 

4 za 


9(J)(x 2 -x^) 


3x 2a D l/l,2^1’ t 2 /n l’ 0;n 2’°^ d? l 


Since, according to (5.37), these are of first order, we need only 
g l/ 2 > E ]_/i 2 " t ° zerot ^ orc ier. From (5.21) and (5.22), 


g l/2 < '^3 ,t 2^ n l’ 0;n 2 , °^ f 1^3’ t 2' > 


(5.43) 
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D l/l,2 (C l’ t 2 /7 b.’ 0; ^2’° ) = «^r^ 1 -v 1 t)+^ 1 t 2 ]) 6C ^i‘^ 1 ) (5 * 44) 


Similarly, in the last term of (5.41), we need only " to zero 'th 

order, which is f^ from (5.21). Putting in (5.41) 


9 gl/1 (t o ,t o /^,0) 




1/12’ 2 l* w/ + v ~ 6 1/1 


at. 


2a 3x 


2a 


o 3 


m 3 v 


2a 


^ a ^ 

3x f l (v 3 ,t 2 )d? 3 f l < ' V 2’ t 2 ) 

2a 


1 3 

+ — 


m 3 v 


2a 




2a 


1 3 ^ ( VVl (t 2' t l )) 


m 3v 


2a 


3x 


2a 


f l (v 2 ,t 2 ) 


(5.45) 


9f x (C,t) 

where we have used — = 0 , so that 

3x 


f^ ) = f x (v,t) (5.46) 

15 

To solve Equation (5.45), we "freeze" f ^ . This assumes that 

gl /1 approaches its t^-* 00 form in a time short compared with that over 

i; 

which f is varying. The solution is 

T; “ A 

Formally, (5.45) can be solved by taking the Fourier transform on x 2 

followed by the standard method of vacation of parameters. It can also 
be solved by using Laplace transform. 
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, 0 ) 


2 2 

2 it k m 


e lk 


3f l (v 2 ) 

3v„ 


2 ik-(v -v )x -ik*v t A . 

e dx + e g 1/1 (k,v 2 ,0/n 1 ,0) 


(5.47) 


where 


-ik"x. 


g l/l(k» v 2» t 2 /n i’ 0) 




( 2tt ) ' 


g l/l < ' ? 2 ,t 2 /n i’ 0)dX 2 (5.48) 


k = k 


(5.49) 


For simplicity, let us omit the term due to initial value. This amounts 
to assuming initially there is no correlation. Equation (5.33) can then 
be evaluated in a straightforward manner, by using (5.47) with X - t -*». 


e < E / 5 ,( 0 ) 
a 1 


V 


= n. 


3<j>(x -x ) 


3 x g l / 1 ^ 2 ,t /^ 1 ’ 0)d ? 2 

la 


n o e 


2m k 


ik ik-(x -X 2 ) Jk i 

2/2 6 g 1 / 1 (^ 2 ,t 1 /n 1 »0)dkdC, 


= n o e 


ik i^-x, - 

dk — — • e 1 (2m) 3 
2m k 


S l / 1^’^2 , t l / \ ,0 )d ^ 
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ik 


ik • x 


= n Q e 


dk 


„ 2 , 2 

2tt k 


1 ( 2tt ) 3 


iff 1 

n 


t =oo 


2 -r£* (v.-v. )t -llc-x.. 

e 2 111 

e 


2 2 

2tt k m 


ik 


8 W 



ik. ( v_-v ) T 
2 1 

e dxdv 


~ n o e 

m 


dk 


k k 
a 


( 2tt ) ° 
(2tt 2 ) 2 


3f l ( ' V 2 ) „ /T >-» -* Jk _> 

— 5 (k- v -k* v 1 )dv 2 

3v 2 


2n 0 e 




(5.50) 


where 


/ k*v 1 

r 

V 

^•v 9 

1 /- 

f i<V 4 

1 

2 

V k / J 

k 

k 


dv. 


(5.51) 


and F' is the derivative of F with respective to its argument. 

We have omitted the initial condition. For the case of non-zero 

initial correlations, there are conditions on the initial correlations. 

18 17 

These conditions have been studied by Sandri and Frieman. ’ 


To find & < E E /?.(0) = n, > , according to (5.34) and (5.37), 

la 16a 11 

we need g and D . to zeroth order, which can be found from (5.21 

Z/_L 

and (5.22). Because of spatial homogeneity, there is no contribution 


due to & 2 / 1 ’ so ' t ^ at 


6 <E la E 16o /5 l <0) = V 
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= ( N-l 


^(x^) 

3<j)(x +v a-x') 

1 1 T) 

3x l<, 

ax i 6 

■ 3 4> Cx 1 -lc 2 ) 

3 4> (3c o-x' ) 

XX X f /• 

9x 

la 

3x lS 1 

3 *(^-* 2 ) 

3<j>(x 1 -x 2 +(v' 1 -v 2 )a) 

3x 

la 

3x lS 


r' 2 -v 2 )6£ 2 dZ' 2 


f 1 (v‘ 2 )d5 2 


n ° 6 J 9 22 

J 2tt k 


ik^ ik* (x^-'x^ ) 8^(jf 1 -x 2 +(v‘ 1 -v 2 )a) 

272" 6 


f 1 (v 2 )dkd? 2 


rr ik (-ik ) -ik*(v-v )a (2ir) ^ 

n e yr-z — e — f.(v )dv dk 

2tt k 2 tt k 


2n„e rr kk„ -ik*(v -vja 


1 2 _ 


f 1 (v^)dv 2 dl<' 


n 1 > da 


2n„e r r r° k..k„ -ik’(v -v 0 )a 


fl(v 2 )dad^dlT 


C r k k , 

2n Q e j I - a ~ 4 - ■ 6 (k- ^-k-T^f^^dv^dk 

, k k. k • v. \ 

n 4 a 6 „ 1 

2n o e ~ ir F — dk 

J k 1 
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Here again we have frozen and taken the limit t -*». 

We are ready to write down Equation (5.30) with all the functionals 
involved explicitly. The following points are observed: 

1. Terms involving <E /£ (0) = "it ><E /if (0) = rf > vanish to the 

XOt 1 X -LpO X X 

order concerned (here second order). For, according to (5.33) and 
(5.37), we need to zeroth order for each average. Because of 

spatial homogeneity, there is no contribution. It follows that 


« E lc.’ E l 6 a /f l <0) ' V ' <E la E l 6 o /f l (0) = V (5 ' 54) 


2. To the order concerned, 


K[E la ,E 16a /C l (0) " n l‘ 1 " <E la E lBa > 


(5.55) 


Since we need & 2/1 and 2/1 (5.34) to zer °th order, which are given 

by (5.21) and (5.22). 


3. tt- < E E > = 0 
3x^ la 16a 


(5.56) 


This is observed in (5,52). 


4 ' 3v <E la E 16a > 

IP 


lP 


3E (x +v a,t +o) 

= <E (x ,t ) I > 

la 1* 1 3v. Q 

lp 

3E (x tvo,t +o) 

* -h/v'i’ 6 L ■ - 

lp 
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3E 

= ° si;- •= VVVW'l 3>t+0) *- 0 < 35T7 E 16» > 

-J-P lp 


!fla p 

0 < 3x 1(J lga 


(5.57) 


because of (5.56). Thus (5.30) can be written as 

3 f, „ a 

•*==■ = J ^7- < E, /MO) = IT, >f, 

3t mdv. la 1 1 1 

1 la 

n 

e 2 3 


3 

da <E E > 1 

( 3v ie 

la 16° 

-oo / 


+ T-^“ T^ - I da < E, E > f 

2 3v n 3v j la 16° 1 

m la 16 J 


- ^ — < E /f (0) = n' >f 
m 3v n la 1 11 

la 


e 2 3 


2 3v. 
m la 


3f, 


da < E n E n „ > 


la 16a 3v 


16 


(5.58) 


8f l 

where we have dropped terms involving pr because of spatial homogeneity. 

dX ! 

Putting (5.50) and (5.53) in (5.58) 


!£i- a 

9t i 9v i 


f V 

r - k 

F i 

k’ v 1 

„ 4 

, 2 V 

,-tkk 

n If • P 

V 

Sf 4 

1 2 

dk — F 

k 

Q. 2 

UJ\ "rT • I 

k 

3 v 

1 m 

k 


m 

k 


1 

u 


(5.59) 


which is the Fokker-Planck equation. 
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5 . 3 The Kinetic Equation 

3 1 3^0 

We now investigate the limit by letting n Q r Q — so that n Q r Q ( — Ty) 

mv o 

= 0(1). (5.35) now gives 

1 : 1 : 1 : e : 1 (5.60) 


We will therefore keep terms up to first order. In the cluster 
expansion, T, and Q are assumed to be of higher order than P, which 
is of first order. As before, we assume spatial homogeneity. We must 
find the various averages in (5.30). 

According to (5.33), we need g-^ to first order for < E la /S l (0) 
rt >. According to (5.34) and (5.60), we need fi- rs t order but 

only the zeroth order of ^ since the term due to ^ is already 
of first order. 

Consider the term involving g 2/ / 1 i n (5.34). We have 


9<t>Cxj 

9x 


-x 2 ) 9<t>(x 1 +v 1 a-x‘ 2 ) 


la 


3x 


16 


g 2/i (f 2 ’ t 1 ; ^3’ t 1 +o/n“i,0)df 2 df ' 


(5.61) 


Using (5.15), (5.16) and (5.17), we find that 


1 

v 2" g 2/l 


(5 2 ,t 2 ;c 3 ,t 3 /c 1 ,t 1 ) 


1 

D(T 1 ",t 1 ) 


{- 


3 

2 TT D(f. ,t.) 
i=l 1 1 



Putting (5.62) in (5.61), 
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2 StK^+v ' a-xp ^ 

n 0 f 9x7 9x~: g 2/l (? 2’ t l’^3 ,t l +a/n l’°' ) d ^ 2 d ?3 

la 13 

34»(5T -3T ) 9^(x +v a-x') _ 

=n 0 1 —^r 5iT~ v ' ”<Vi*3-V» dt 2 dq 


2 3 (iTj^ -"x 2 ) a^x^+v^-xp 

n n J jw 


f 1^3 t l +a)g l/1^2 ,t l / ^3' t+a)d ^2 d ^3 


(5.63) 


Here we have made use of the fact that integrals of the type 


3<t>(x n-X _) 

J jj- WV dt 2 = 0 


(5.64) 


and may be omitted because of spatial homogeneity. 

Equation (5 . 63) together with the fact that we need D . to zeroth 

Z 9 Z / -L 

order reveals that 


' E l ( , I Wh< 0) = < E la E lga > 


(5.65) 


and this average is a functional of g-^/i* Thus it will be sufficient to 
determine g^_/± fi rs ' t order. Also, because of spatial homogeneity, 
terms involving products of the average of E^ (conditional or not), may 
be omitted to the order concerned. Hence Equation (5 . 30) may be rewritten 


3f 2 0 

= — -r < E /c7(0) = ‘n. > f, - — V- — j do(-o) 

9t, m 9 v , la 1 112 9v, J 

1 la m la 


< 9^ E 13c > f l 
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e 3 3 


2 3v 3 v 
m la 13 -■ 


I da < E 1 E 1R„ > f 1 + -£r- I < E 1 r* E 1 Rm >f 


la 13a 1 2 3v, 3x., 0 

m la 13-o 


la 13a 1 

(5.66) 


where 


3<j>(x -x ) — _ 

6 < E la /5 1 (0) = ^1 > = n 0 / — if- gl nCto.Vn-, ,0) de„ (5.67) 


3x, to l/l 2 ’2 1 : 

la 


according to (5.33) and 


e < E, E n > 
la 13o 


n n I 


d<p(x -xT 0 ) 3f(x.+v a-x') 

. f l U 2'V 


3x 


la 


13 


6(v’-^ 2 ) d? 2 dq 


34> (3< -x , ' ) 3i(i(? +v 3-x') ^ _ 

+ n 0 J — — g— f l(^3 ,t 1 +CT ) gl/l U 2 t l /5 3’ t l +0) dC 2 d ^3 

la 13 

(5.68) 


according to (5.34), (5.63) and (5.65). Here we have substituted 


V f i (? 2 ,t l ) 


6 ( x ^- x 2 - v ‘ 2 o ) 


6(v’-^ 2 ) = 


D 2,2/l (? 2 ,t l’ 5 2 ,t l ta/n l’° ) 


to zeroth order. 

To find g ^ to first order, let us rewrite Equation (5 . 32) , omitting 
products of averages because of spatial homogeneity. 

9g l/l ( ' ? l ,t l / ^2 ,t 2 ) 3g l/l 

1 + v — 

3t, la 3x, 

1 la 



71 


m 9v, < E la /5 1 (0) V ^2 (t 2 ) "2 > g l/l 


2.0 3 E 

“2 9v~ $ da(-o) < = ’V 5 2 (t 2 ) ? 2 > g l/l 

m la 


3 x lg 13a^r 


e 2 3 


2 9 v 9 v 
m la 13 


/ d0 ' E iaW f i <0) -- E • MV =t, > g 


1 ’ 2 2 ' 


’ 1/1 


e ^_ _9 9 _ 

m 2 9v la 9K L3 


| da(-a) < E Z /KAO) = n , , t,(t ) = t 2 > g 


J la’“16a' ’’I"" "1’ ’’ 2 ' “ 2 ' ^2 ' & 1/1 

(5.69) 


Now 


e < E la^l^°^ n l’ ^2^ t 2' ) ^2 > 


N 9<j>(x -x.(t )) ^ 

= < E 31 T - = n i ’^ 2 ( t 2 ) ? 2 > 

j =2 la 


_ f !! a ^ Oyx .) g 1/2 

j=3 9x la 


9x, V (5j, t 1 /n 1 ,0;C 2 , t 2 ) <U.. 


S ^ Cx ^- x ’ ) 

jf 9 Z 
la 


D 2/l,2^2’ t l /n l’°’^2 ,t 2 ) d ^2 


9 <))( x 1 - x 3 ) 

= n 0 J” 9^T B l/2' s 3 , -l' "i’“’-2’^2' "“’3 

la 


g n /0 (T q ,t /tu ,0;£ ,t ) d£. 


94>(x -x’ ) 


+ J — - 3 ~ -— fi(x’-C^ 2 +v 2 (t 1 -t 2 )]) 6 (vVv 2 ) dC’ 


la 


(5.70) 


where we have substituted the zeroth order of ^ 2/1 2 s ^ nce "term is 

of first order already. 

Again we use the cluster expansion to find g-^- Using (5.16), 
(5.17) and (5.64) , we find 
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3(|>(x -x ) ^ ^ 

n 0 i" 3x. g l/2^3’ t l /n l.’ 0: ^2 ,t 2 > ' d ^3 

la 


= n. 


3<(,(x 1 -x 3 ) 


— g, /-i Uo>t,/n, ,0) d£. 


3x, to l/l 3 ’ 1 1 : 

la 


+ n. 


a^C^-Xg) 


3x 


la 


g l/1^3 ,t l^ ri 2 :t 2' > d ^3 


so that (5.70) becomes 


(5.71) 


6 < E la /? 1 (0) = ^i’ 5 2 (t 2 } = 5 2 > 


= n„ 


+ n. 


g l/1^3 ,t: i /r| l ,0) d? 3 


3<f)(x 1 -x 3 ) 

9x 

la 

8<()(x 1 -x 3 ) a 

3x~ g l/l (C 3 ,t l /n 2 ,t 2 ) d£ 3 

la 


3<Kx 1 -x 2 -v 2 (t 1 -t 2 ) ) 


3x 


(5.72) 


la 


Next, we want to show 


< E la E iea /5 l (0) = V ? 2 (t 2 } = ? 2 > 


E la E lBo 


(5.73) 


For 


e2 < E la E 16o / h^ 0) = VVV = 52' 


N 3(j)(x, -jf. ( t, ) ) 9<()(x 1 +v, a-'x, (t +a) ) 

< E — i-^L — /«!«» 


]<k 

j ,k>2 


la 


IB 


"^l’^2 (t 2 ) = ^2 > 
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N 3<j>(x -x. (t ) ) 3<j>Cx +v“ .o-x (t. +a) ) 

+ < E — k - 1 S- x ------- — /r x (o) = vt 2 (t 2 ) =t 2 > 


+ < z 


j =k> 2 ““la ”"'13 

N 8 1(> 0^ ( t - l ) ) 3<t>(^ 1 +v > 1 a-'x (tj^+o) ) 


k >2 ^ X la 


3x 


13 


•/C 1 ( o) = n 1 ,5 2 (t 2 ) - S 2 > 


N 9<f>(x -x. (t )) 3<t>(x +v. a-'xit .+a)) , . 

+ < . s 9 — k ; 2 — ^i (0) = V 2 (t 2 ) = « 2 > 

] > 2 la 13 


3<|)(x -,-x* ‘(t )) 9<j>Cx tv o-$ (t +a) ) ^ ^ 

+ < — k 1 ^ kr- 2 - — ^1 (0) = WV = ^2 > 


la 


16 


(5.74) 


The last term is of second order. The third and fourth terms are 
proportional to N and are of first order, so that in evaluation, we need 
only density functions to zeroth order. To this order the average on 

-_y t 

x,(t +a) or on x.(t ) (j,k > 2) has no contribution because of (5.64). 

K _L 2 -L 

In the second term, the density function needed is of zeroth order and 
in fact this term is the same as the first term of (5.64). For the 
first term, we have 


N 3<j>(x -x.(t .)) 3<J)Cx +v o-x, (t +o) ) 
„ r 1 3 1 11 k 1 

.,r 3 x 

3 ?k la 

j ,k>2 


3x 


16 


•/£ 1 (0) 


n l , ^2^ t 2^ ? 2 > 


3<f>(x -x ) 3<{)0< +v a-x ) a 

- n„ J ^ g 2 ^ 2 (^ 3 ,t 1 ;C 4 ,t 1 +a/n 1 ,0-,C 2 ;t 2 ) d^,^dE, L 


3x 


3x 


la ”“13 

2 9<()(x 1 -x 3 ) 3 <|> x ct - x 4 ) 

_ n r 


0 


j 


3x 


la 


3x 


13 


f i ( 5 4» t 1 +a) g l/l ( 5 3’ t l /5 4’ t l +a) 


d « 3 < 


(5.75) 
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by using the cluster expansion and the spatial homogeneous property 
(i.e. (5.64)). Essentially the reasoning is as follows: 


D 


v 2 g 2/2 


1 , 2 , 3,4 

° 1,2 


Because of (5.64), we may restrict D to its first order component, 

^ 1 , 2 , 3 ,4 

which is £ D.D.D , and D to its zeroth order component which 

• • i ^ 1 K » Jw -L * z. 

K] ,k<£. J ’ 

is D^D^. Again because of (5.64), the only surviving term is when 
k = 3, l = 4. 

Thus (5.73) is valid. Since (5.72) and (5.74) are all of first 
order, we need to zeroth order on the right hand side of (5.69), 

which is f . Substituting into (5.69), 


9g l/l ( ^l’ t l /? 2 ,t 2 ) 3g l/l 


at. 


+ V 


la 9x 


la 


n 0 9 . a4>(x-,-x Q ) 

m 9v n 
la 


9x g L/l ( ^3’ t l / ^2 ,t 2 ) d5 3 f l (5 l 9t l ; 

la 




m 9v, J 9x 6 1/1 

la la 


E (f 3 ,ti / Tti.°) df 3 f 1 (f 1 ,t 1 ) 


i a ^(yyVytp) 

m 9v, 9x, 

la la 




e ! 3 1 3E lc 

J do (_ o) < 

lp 


2 9v 
m la - 1 
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e 1 2 * 3 3 f 

m 2 3v la 9V 16 i ° 


E la E lga > f l (C l ,t l ) 


2 3v 3x * d ° ( o) * E la E lga > f l (C l’ t l ) 

m la lg 


Finally, using (5.66) 


9g l/l^l ,t l /? 2 ,t 2^ 


la 3x, 


’o 3 , 8 * (x l- x 3 ) 


u r 

m 9v ' 


3x g l/l < ' S 3 ,t l /? 2 :,t 2' ) d? 3 f l < ' ? l ,t l ) 
la 


+ 1 9 4 > ( x 1 -x 2 -v 2 (t 1 -t 2 )) 

m 3v, 3x_, 

la la 


— f i (5 1 ,t 1 ) + 


9 f l ( 5 l’ t l ) 


Since we are interested in the functionals of g such as < E (0) 

— f- -L / -L _LCX -L 

df i 

a > and < E, E, „ > . we may omit the term 7 -— in (5.77) which has no 

1 la 16a ' 3t 

contribution because of spatial homogeneity. (This amounts to writing 
= f^ + 6 g and consider only 6 g). Furthermore, we ,: freeze ; f^ as 
before, in finding these functionals, so that 


9 g l/l (g: i ! ’ t 1 ^ 2 ;t 2 ? 9 g l/l 

3t, + V la 3x 

1 la 


n n 3 3<f>Cx -x ) , ^ ^ 

in "377“ ^ g l/l (C 3’ t l /C 2 ,t 2 ) d? 3 f l (v l } 


1 3 9< * >( V" 2 -"2 (t l- t 2 )) 

+ in 3v~ 337 W 

la la 
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By symmetry , 


3g l/l ( ' ? 2’ t 2 / ^l’ t l^ 3g l/l 

8t 2 + V 2a 8x 2a 


^0 _J r T 

m 9v_ 9x 


34>(x 0 -x q ) 


„ g l/l (5 3 ,t 2 //5 l ,t l ) d? 3 f l ( ^2 ) 


1 3 3^x 2 -Xi- v i (t 2-ti )) 

+ m 9v 9x~ 1 V 2 

2a 2a 


(5.79) 


so that 


3g l/l (g 2’ t 2 /n l’° ) 9g l/l 

9t 2 V 2a 9x 2a 


n_ . 94>(x -x ) ^ . 

“ lw7~ i" ~9^r g l/l (5 3’ t 2^1’ 0) d ^3 f i (v 2 ) 


i 3 


+ m 9v, 


— f^Vj) 


(5.80) 


Equation (5.80) may be solved by using Laplace and Fourier transforms. 


g l/l (k 2 ,V 2 ,p) 


00 "Pto i o -v _v 

1 e dt (j-) / e gi/i ( ?2’ t 2 /n l’° ) 

0 


(5.81) 


and define the operator 


V(2) = 


,3 9f , ( v 0 ) ik o <t>(k 0 ) 


, sG Oi (V n ; -LJX Y \ / 9 

L 2 i. ) ...qo. - p -l J. — 2 r — --r-w— ] 

m 53 E ( k97d „ ) J P 2 + lk 2 *V 2 


(5.82) 
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2 2 
2tt k‘ 


(5.83) 


is the transform of <(>(x), and 


€(k 2 ,p 2 ) = 1 E 


( 2tt ) n Q <j)(k 2 ) 


9f (v ) dv 

J ^ ■ ■ ■ ! ^ - (5.84) 

8v 2 P 2 + ^2^2 


Sl/fVVM * V(2) S 1 


(5.85) 


1 8f l (V 2_l 

m 


ik 2 (|)(k 2 ) -i1< -(x^-v!^) 

— . a • - e 


P 2 + 1 ^2' V 1 


(5.86) 


is the transform of — v - • - 

m 9 v 2 


1 3 H(k 2 -x 1 -^ 1 (t 2 -t 1 )) 


f i<n> 


Here we have ignored the initial condition. We find that 


I S in (k 2 ,%,P 2 > iv 2 




e(]< 0 ,p 0 ) 1 p 2 + i^ 2 ''^2 m 8V 


i^ 2 4>(k 2 ) -iT< • (3^-v.^) 

7—k — =* — e 


p 2 + ^1 


Taking t 2 = t -* 00 , we may evaluate the residue to obtain 


f g (k ~v ,t = t -*■ °° ) dv„ 
J S l/1^ 2’ 2’ 2 1 2 
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- ik -x 

p 1 , .x , , 

-r* -T* -o (1 - e(k.,-ik -v )) 

ed^.-i*^) (27i) 3 n 2 2 1 


0 


(5.88) 


From (5.67) 


e < E / ?1 (0) = n, > 


a*(x 1 -x 2 ) 

n 0 t 3jT e l/l 

la 


g i(f 2 , t i / ni, 0) df ; 


ik-x ^ 

n Q / dk ik^4>(k) e (2tt) / k > v 2 ;t l^ dv 2 


= J dk ik a *(k) 


(1 - e(k,-ik*v 1 )) 


/ dk y(k) Im el¥ -_-^-y 
( 2tt ) 3 n e 2 

k e 2 -rrrr-^' ( - ^ } 

r a 2it mk k 

i dk— o 


where 


2iT / k^ | e(k ,-ik>v 1 ) | 2 


2n 0 e 

m 


^ -A 

k-v 

k F ' (— r— — ) 
r dk a k 

J ~TT . , jk -x -x . | 2 

k j e(k ,-ik • v^) | 
(5.89) 


t- 


■A -r»> Jk _» 

v k-v k-v 

F(— ) ° / f(v 2 ) 6(— - ^ — ) as before. 


We now turn to the correlation term. According to (5.68) we need 
f 1 (q,t 1 + c) g 1/1 (t 2 ,t 1 /f’,t 1 + o) = V 2 D(f 2 ,t i; t;,t to). In other 
words, we need a two-time two-body distribution. Multiplying (5.79) by 


f 1 (? 1 ,t 1 ), 
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3 V D(5 2 ,t 2 ;C 1 ,t 1 ) 


3t, 


+ v 


aV 2 D(t 2 ,t 2 ;'f 1 ,t 1 ) 


2a 


3x 


2a 


JO J r 

m 3v J 


S*<x 2 -x 3 ) 


m 3v 


2a 


3x 


2a 


V D ^3’ t 2’^i’ t i^ d ^3 f l^2^ 


3<f)(x -1c -v“ (t -t ) ) 

,13 21121 ,-*■ s,- \ 

+ m 3^ 9^7 f l v 2 f l V 1 


2a 


la 


(5.90) 


, s, .y n y 

Here we have omitted the fjC ? 2 ,t 2 ) • f^ ^ jt.^) component of V 1(5 ,t ; 

^l ,t: i 3 us ’ t as dn (5.79) since it has no contribution to the correlation. 
Equation (5.90) is to be solved subject to the initial condition 

V 2 D(t 2 ,t 1 ;f 1 ,t 1 ) at t 2 = t ± . 

2 jk ^ 

To find V D( C 2 »t^ ,t ) , we add to (5.90) its symmetrical counter 
part by interchanging the subscripts 1 and 2, setting t 2 = t^. The 
result is 


3V 2 D(f 2 ,t 1 ;f 1 ,t 1 ) 2 3V 2 D(f 2 ,t 1 ;? 1 ,t 1 ) 

3t + V ja 3x . 

1 :=i :a 


n n a 3<t>(x -x . ‘) 

= — — / — v D(C 3 , t i;? 2 , t i ) d5 3 f l ( ^l } 


m 3v n 1 3x 3 

la la 


n n a 3<()Cx,-x ) ^ . 

+ I — V 2 D(5 Q ,t 1 ;C 1 ,t 1 ) d5 0 fjvj 


m 3v J 3x 

2a 2a 


’3 1 1 ’ 1 3 1 2‘ 


3*(Xi-x 2 ) s SMXj-xp 
m C 3v, 9x, 3v„ 8x„ ] 


(5.91) 


la la 


2a 2a 
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2 ^ 

For t sufficiently large so that V £^,t^) r e l axes to its 

asymptotic form, (5.91) may be solved omitting the initial condition. 

^ a 

This is then used in (5.90) to solve for D(C 2 5 t 1 + which now 

2 

depends only on x. Substituting this into (5.68), e < E la E 1 g 0 > can 

21 

be evaluated. This has been done by Rostoker. The result is 


< E la E 16o 


c da) ima 

J 2tT B 


a i . a) , 

, dk ik 'V , u ,2 2 k o k B 2. F< - k> 

J e ( 4tt ) n„e 


( 2tt )" 


0^ , 4 k | . s 2 

k e ( k . iu) ) 


Integrating over c, 


(5.92) 


/ d ° " E la E 16a " 


= J 


dk 


( 2tt )' 


k k p ( - — ) 

r da> Aa. ,2 2 a 6 2 tt k' 

J — it 6(a) + k-v, ) ( 4 tt ) n n e 


2 it 


■°' k 4 k le 


0<,ia)) 


k k, 
a_ 

C 

k" 


H 

e( -t A) 


o 2 r a t a 

= 2 "o e I dk T5 , A - - ,,2 


e(k»-ik*v ) 


(5.93) 


Since < E E . > does not depend on 1, Equation (5.57) is valid, 

la 10a 1 

Furthermore, the last term in (5.66) vanishes. Putting (5.89) and 


(5.93) in (5.66), we obtain 

4 


k • v. 


3f , 


3 2 V 


3t 3v 2 

1 la m 


k F 

f j, a k j- 

I dk ~u f 


k | e(k, -iT<- v^) 


2 1 
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3 ( 3 


2n„e 


k k. 


3v 3v 2 
la 18 m 


° w [ 'of'B 
dk 


k-v 

r<- k -A) 


5 , o. _v i o ^ ^ l 

k | e(k., -ik.v^) j 


3v 3v 2 J 
la 18 m 


"t .~ s " 

4 k V ] 

2n e k k F(— - — ) 

0 r d J -JLl k 


k 5 I e(k ,-ik-v^) I 2 1 


3v„ 


2n n e v * 

0 f jT k 

o j dk 


f k-t 

F( ^ 


v m 


4 ,-jw a v 9 1 1 
k | e(k,-ik ■ v ) | 


2n e 

+ “I— J dk 

m 


v“ ^ 
k-v 

■*kk k ' 


3f 


k 5 | efk.-ik-v^) | 2 3v l 


1 


(5.94) 


k ' V J 

Letting u^ = — 


, (5.94) can be rewritten as 


3f ] 

3t, 


3 

3ir, 


2n 0 e 


m 


, -*■ 

/ dk 


k 


[fn (V n ) 


3F( V 


k 4 | e(k, -ilT-'v ) | 2 1 1 9u l 


3f , ( v. ) 

-F(u ) X ■ - ■ ] 
1 9u, 


(5.95) 


Equation (5.95) is the kinetic equation. 


19 


5.4 Discussion 

We have seen in the previous sections, decoupling the set of the 
generalized stochastic equations by using the cluster expansion and 
freezing the lowest distribution, result in the Fokker-Planck or kinetic 
equation. We have assumed that the electric field experienced by a 
particle is of small order. This of course cannot be proved, and is 
simply false when particles are close. There is, however, another 
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problem. Actually, the force experienced by a particle consists of 
forces due to a large number of other particles, and we have only assumed 
each individual force is small. In other words the electric field 
seen by the jth electron is given by 


^ N ^ 

E. = £ E. . 

3 i*j 13 


where 


a 
E. . 


0(e) 


We wish to indicate it is reasonable to take 


Ej = 0(e). 


With the Fokker-Planck case, where n.r„ = 0(1), E. = 0(e) since the 

0 0 j 

3 

summation actually only extends to N ^ n 0 r 0‘ With ’ t ^ e kinetic case, 
3 1 

however, n^r^ □ 0(— ) . Let us write 


N F. 


ID 




1 -V 

where F. . = 0(1), rr = 0(e). To zeroth order the F..'s are identically 
i] N i] 


distributed and independent, since the trajectories of the particles 
are independent and have identical distribution. Thus E^ -*■ < F^ > = 0, 
as N ->-oo, which indicates E_. is small. 

When the cluster expansion was introduced, we required T, Q of 
higher order than P. Actually this is not necessary for the case 
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n^r^ = 0(1). Also in this case, the omission of the initial correlation 

is reasonable in most cases. In fact the contribution of the initial 

correlation as given by (5.47) can be written as g^O^-v,^ ,v* 2 >°/ 

"x-^-v^t "v^,0) with ± 2 n -*■ «>. This will be small if the difference 

between the spatial arguments is large, which is the case with the 

exception • This is expressed by Sandri as '‘the principle of 

18 

absence of parallel motion.' 1 

For the kinetic case, the major problem is in the evaluation of 
the correlation of the electric field since this involved the asymptotic 
distribution of two particles. The following, however, is observed. 

If the equilibrium distribution is used, the correlation is found to be 
formally the same as given by (5.92), 


< 


E la E 16a 


> 


r 

j 


da) 

2tt 


iaicr 

e 


d" ik ' V l° (4tt) 2 Vb 2tt 

( 2tt ) 3 k 4 k 


F(-£) 


.a 2 

e(k ,iai) 


where 

F(u) = / f (v) 6(u - ~~ )dv' 

J m k 

and f is the Maxwellian distribution. If the distribution function 
m 

f^ differs from f only in first order, then (5.92) may be expected 

since only the zeroth order of F is needed. 

Another way to obtain the fluctuation of the electric field is to 

-4 22 

make use of the concept of the dielectric constant e(k,ioi) or the 
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concept of dressed particles. While these techniques yield almost 

immediate result, a mathematical justification, however, involves 

23 

considerable difficulty. 

a 2 

The term < E.^/5^(0) = > gives the polarization effect, also 

9 

referred to as the field due to quasi-particles by Rostoker. Note 

that in the kinetic case, this has an extra factor w 

|e(k,-ik^)| 2 

compared to that in the Fokker-Planck case and in fact this factor is 

the only difference between the kinetic and the Fokker-Planck equation. 

2 9 

This force has been obtained by Gasiorowicz, et al. and Rostoker by 

using the linearized Vlasov equation. The use of the Vlasov equation is 

essentially justified, according to (5.80). However, in Rostoker : s 

calculation, the Vlasov equation is solved with either an external 

driving force or a particular initial distribution, so that there 

seems to be some arbitrariness. When an external force is assumed, it 

is taken to be that due to an electron moving in its rectilinear orbit. 

-i£. v^t 

The resulting polarization force has an extra factor e , which is 

not explained. Such a factor does not appear in our calculation and is 

due to the condition 5,(0) - n, in < E, /5,(0) = n, >• Here n, = 

1 l la l l x 

(x^-v^tjV^). If the condition had been £^(0) = 5^ = Cx^'v’^), t ' rie 
factor e 1 will appear. Another point to be noted is that while 
we obtain the polarization effect, no test particle is introduced 
artificially . 

The Fokker-Planck and the kinetic equations (5.59) and (5.95) have 

17 19 

been discussed in the literature. ’ It is to be emphasized that we 
do not obtain these equations from the BBGKY hierarchy, although the 
mathematics turn out to be similar. This similarity has deeper signifi- 
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cance than it appears, as we will show in the next chapter, and justifies 
to some extent the hybrid use of the conventional Fokker-Planck 
equation and the BBGKY theory. The generalized stochastic equations, 
of course, are applicable to very general situations. 
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VI. FINAL REMARKS AND CONCLUSION 

In the previous chapters we have shown how a set of generalized 
stochastic (or generalized Fokker-Planck) equations can be used to 
describe the evolution of a system. When applied to the Coulomb 
potential problem, a strong similarity between the set of generalized 
stochastic equations and the BBGKY hierarchy is observed. In fact, the 
BBGKY hierarchy may be regarded as a special mode of generalized 
stochastic equations. To see this let us take the limit t -*■ 0 in 
Equation (2 . 3 1). All terms involving integration vanish, leaving 


9u) 

9t 


t = 0 


9 _ 

9x 


< e F(x,0)/X(0) = x > to 


t = 0 


Since we have chosen 0 arbitrarily as the initial time, in general, we 
have 


|qr = - < e F(x,t)/X(t) = x > w 

9t 9x 

This equation is exact and a set of generalized stochastic equations 
obtained in this manner is equivalent to the BBGKY hierarchy/' Thus, 
in this sense, the BBGKY hierarchy may be regarded as a special mode 
of generalized stochastic equations. Note that while it is not possible 
to ignore the condition X(t) = x in < F(x,t)/X(t) = x > for all t, 
since X(t) depends on F, it is possible to ignore the condition X(0) □ 


See appendix 
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x in < F(x,t)/X(0) = x > provided there is no correlation between 
F(x ,t ) and X( 0 ) . 

Usually the Fokker-Planck equation is derived by using a Markov 

2 24 

model or using the Lagrange expansion. ’ The coefficients are 
averages of the force and no explicit condition is imposed. We have 
seen in Chapter V that while we can drop the condition on the correlation 
of the electric field, the exact condition is important in calculating 
the polarization effect. In fact for a spatially homogeneous plasma, 
the unconditional average of the electric field will be zero. 

Our results may be summarized as follows. We modify Stratonovich ' s 
results to obtain a set of generalized stochastic equations which 
describe the higher distributions as well as the lowest forming a chain 
similar to the BBGKY hierarchy. It is possible to obtain the kinetic 
equation by decoupling this chain with the cluster expansion. This set 
of generalized stochastic equations has the following characteristics 
which distinguishes them from the conventional Fokker-Planck equation. 
Firstly, the coefficients involve conditional averages of the force 
taking into account the correlation of the force and initial position. 
Moreover, the equations serve to describe higher distributions as well 
as the lowest distribution, making it possible to determine the 
coefficients without introducing certain models or appealing to other 
theories. There is a similarity between the generalized stochastic 
equations and the BBGKY hierarchy. In fact the BBGKY hierarchy may 

be regarded as a special mode of the generalized stochastic equations. 

1 2 

For a plasma Sturrock, Gasiorowicz, et.al. obtained essentially 
the same Fokker-Planck equation. To determine the coefficients, 
Gasiorowicz, et.al. use the Vlasov equation and a Holtsmark distribution. 
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This may be inconsistent although we find the use of Vlasov equation is 

essentially justified. Also we find that an artificially introduced 

test particle is unnecessary. In an entirely opposite direction, a 

kinetic equation on the basis of the BBGKY theory may be manipulated 

9 

into the form of a Fokker-Planck equation as done by Rostoker and 
Tchen.^ Our approach is to regard the set of generalized stochastic 
equations as basic from which the coefficients in the equation for the 
lowest distribution are determined resulting in the kinetic equation. 

In this case either approach involves similar mathematics. However, 
the generalized stochastic equations are applicable for general 
situations. We have considered several situations where the random 

force is assumed to be specified, with some results comparable to those 

1 7 

of Sturrock and Puri. 

We conclude that the generalized stochastic equations give a 
general formulation applicable to self-consistent systems as well as to 
systems with specified random forces. 
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APPENDIX 


We give a heuristic argument to obtain the equation 


9 (u 9 

9t 9x 


< F(x,t)/X(t) = x > 


(A.l) 


The notations of Chapter III will be adopted. 
Consider the system, 


x = F(x,t) 


(A. 2) 


Let w(x,t) be the probability density of X(t). Firstly, we assume 
F(x,t) is a given function. We have 


w(x 0 ,0) = iu(iHt,t 0 ,x 0 ),t)j 


(A. 3) 


where J is defined by (3.2). Differentiating (A. 3) with respect to 
t, we obtain 


coOJj ,t)J = 0 


or 


Tr 9m(^ ,t ) , ^ \ ,t ) -, , ,, ^ dJ 

J[ at + F j ( ^’ t) Tx. ] + dt 

: 


= 0 


(A. 4) 


Now 


dJ T 9F. 

at = J -J- 

9x . 
J 


(A. 5) 
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along the solution. Putting in (A. 4), cancelling J, and replacing 4* by 
x, we have 


9w(x,t) 

3t 


3x , 


Fj(x,t)o)(x,t) = 0 


(A. 6) 


If the force F is random, we can write Equation (A . 6) for each 
sample F; here w(x,t) should be replaced by u)(x,t/F), the conditional 
density of X(t) given F. Thus, 


9oj(x,t/F) 

at 


F.(x,t)u>(x,t/F) = 0 
-n i 


Averaging over the function space of F, 


(A. 7) 


+ gj- J F.(x,t)oj(x,t/F)dP(F) = 0 (A. 8) 


where P(F) is the distribution of the force F. Formally we take 


u>(x ,t/F )dP( F ) = w(x,t)dP(F/X(t) = x) 


(A. 9) 


according to the usual rules on conditional probability. Here 
P(F/X(t) □ x) is the distribution of F given X(t) = x. (A. 8) then 
yields 


/ F (x,t)dP(F/X(t) = x) m(x,t) = 0 

j 3 
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5ai(x,t)- + 3 — < F.(x,t)/X(t) = x > u(x,t) = 0 

dt dX . I 

] 


which is (A.l) 


Let us apply the result to (5.1) with j □ 1. (A. 6) gives 


H 1 ' 51 ’ 41 + V t)/? (t) = V V*» = 'l * f l 


+ -Z~- < ^ - 8 ^i yx (t) = x , V (t) = \ > f x = 0 

3 V.. ~ ~ 1 

1 \ 3x l I 

1 (A. 10) 


< v^Ct) /f ± (t) = \ t v^Ct) = 


< ~ (- 9 ^l|/'^ 1 (t) =‘x 1 , V-^t) = v JL > 


< E - (x i X j (t)) /x (t) = x , V t .(t) = v 

jtl m 3x x X 1 11 1 


E .If v 2 f2(?1 ’ t} d5j 

j+1 m 3x 7 T 

V f l (5 l» t} 


1 N-l r 9<() ■* x ^ /■? -f . \ 1 

TT—J (^-i ) -^9 ' ^ -i 5 ^-^9* v 

m V J ^ 1 z z l z l - 


fjU.t) 


Putting in (A. 10), 
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+ . 50 jl 

9t m 

9x 3v^ 


I s^_ (? 1 - «,) f (t r tut 

9x^ 


(A. 11) 


which is the first equation in the BBGKY hierarchy. The entire hierarchy 
including the Liouville equation can be obtained in a similar manner. 
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Vlasov equation in finding the polarization effect] is esjsentiajlly justified. Ttje 
polarization effect is seen to be the consequence jof the: condition ilmposed on the 
average in evaluating the coefficient. Unlike the test particle thejory, however, 
no test particle is artificially introduced. Finally, a similarity jis observed Ibetwe^n 
the set of generalized stochastic equations and thle BBGKY hierarchy.! Closer 
examination reveals that the BBGKY hierarchy may be regarded, jin a dense, 'as a ^pecia 
mode of generalized stochastic equations. 1 j 


Security Cln-ctui -.i-i n 










